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Brane-world cosmology is motivated by recent developments in string/M-theory and offers a 
new perspective on the hierarchy problem. In the brane-world scenario, our Universe is a four- 
dimensional subspace or brane embedded in a higher-dimensional bulk spacetime. Ordinary matter 
fields are confined to the brane while the gravitational field can also propagate in the bulk, and 
it is not necessary for the extra dimensions to be small, or even compact, leading to modifications 
of Einstein's theory of general relativity at high energies. In particular, the Randall-Sundrum-type 
models are relatively simple phenomenological models which capture some of the essential features 
of the dimensional reduction of eleven-dimensional supergravity introduced by Hofava and Witten. 
These curved (or warped) models are self-consistent and simple and allow for an investigation of 
the essential non-linear gravitational dynamics. The governing field equations induced on the brane 
differ from the general relativistic equations in that there are nonlocal effects from the free grav- 
itational field in the bulk, transmitted via the projection of the bulk Weyl tensor, and the local 
quadratic energy-momentum corrections, which are significant in the high-energy regime close to 
the initial singularity. In this review we investigate the dynamics of the five-dimensional warped 
Randall-Sundrum brane-worlds and their generalizations, with particular emphasis on whether the 
currently observed high degree of homogeneity and isotropy can be explained. In particular, we dis- 
cuss the asymptotic dynamical evolution of spatially homogeneous brane-world cosmological models 
containing both a perfect fluid and a scalar field close to the initial singularity. Using dynamical sys- 
tems techniques it is found that, for models with a physically relevant equation of state, an isotropic 
singularity is a past-attractor in all orthogonal spatially homogeneous models (including Bianchi 
type IX models). In addition, we describe the dynamics in a class of inhomogeneous brane-world 
models, and show that these models also have an isotropic initial singularity. These results provide 
support for the conjecture that typically the initial cosmological singularity is isotropic in brane- 
world cosmology. Consequently we argue that, unlike the situation in general relativity, brane-world 
cosmological models may offer a plausible solution to the initial conditions problem in cosmology. 



1 



I. INTRODUCTION 



Cosmological models in which our Universe is a four-dimensional brane embedded in a higher-dimensional spacetime 
are of current interest. In the brane-world scenario, ordinary matter fields are confined to the brane while the 
gravitational field can also propagate in the extra dimensions (i.e., in the 'bulk') [4-6,192,204]. In this paradigm 
it is not necessary for the extra dimensions to be small or compact, which differs from the standard Kaluza-Klcin 
approach, and Einstein's theory of general relativity (GR) must be modified at high energies (i.e., at early times). At 
low energies gravity is also localized at the brane (even when the extra dimensions are not small) [7,41,193,194]. 

The five 10-dimcnsional (10D) superstring theories and the 11D supergravity theory are believed to arise as different 
limits of a single theory, known as M-theory [123,191,211]. The 11th dimension in M-theory is related to the string 
coupling strength, and at low energies M-theory is approximated by 11D supergravity. In Horava-Witten theory [115], 
gauge fields of the standard model are confined to two 10-branes located at the end points of an S 1 / Z2 orbifold and 
the 6 extra dimensions on the branes are compactified on a very small scale effectly resulting in a number of 5D 
"moduli" fields. A 5D realization of the Horava-Witten theory and the corresponding brane-world cosmology is given 
in [157-159], with an extra dimension that can be large relative to the fundamental scale, providing a basis for the 5D 
Randall- Sundrum (RS) models [193,194,92]. An important consequence of extra dimensions is that the fundamental 
scale is M4 + d, where d is the number of extra dimensions, rather than the 4D Planck scale M p = M4. If the extra- 
dimensional volume is significantly above the Planck-scale, then the true fundamental scale M^+d can be much less 
than the effective scale M p ~ 10 19 GeV and the weakness of gravity in 4D can be understood in terms of the fact 
that it leaks into the extra curved or "warped" dimensions. 

As in the Horava-Witten theory the RS branes are Z 2 -symmetric and have a self-gravitating tension, which counter- 
acts the influence of the negative bulk cosmological constant on the brane. The RS brane-worlds and their generaliza- 
tions provide phenomcnological models which capture some of the features of M-theory, they provide a self-consistent 
and simple 5D realization of the Horava-Witten supergravity solutions [115], and they allow for an investigation of 
the essential non-linear gravitational dynamics in the high-energy regime close to the initial singularity when moduli 
effects from the extra dimensions can be neglected. RS brane-world models have in common a five dimensional space- 
time (bulk), governed by the Einstein equations with a cosmological constant, in which gravity is localized at the brane 
due to the curvature of the bulk (i.e., "warped compactification" ) . A negative bulk cosmological constant prevents 
gravity from leaking into the extra dimension at low energies. In the RS2 model there are two ^-symmetric branes, 
which have equal and opposite tensions. Standard model fields are confined on the negative tension brane, while the 
positive-tension brane has fundamental scale M and is "hidden" . In the RSI model there is only one, positive tension, 
brane which can be obtained in the limit as the negative tension brane goes to infinity. We concentrate here mainly 
on RSI branes in which the negative cosmological constant is offset by the positive brane tension; the RS2 brane 
models introduce the added problem of radion stabilization and complications arising from negative tension. 

It has recently become important to test the astrophysical and cosmological implications of these higher dimensional 
theories derived from string theory. In particular, can these cosmological models explain the high degree of currently 
observed homogeneity and isotropy? Cosmological observations indicate that we live in a Universe which is remarkably 
uniform on very large scales. However, the spatial homogeneity and isotropy of the Universe is difficult to explain 
within the standard general relativistic framework since, in the presence of matter, the class of solutions to the Einstein 
equations which evolve towards a Friedmann universe is essentially a set of measure zero [63]. In the inflationary 
scenario, we live in an isotropic region of a potentially highly irregular universe as the result of an accelerated expansion 
phase in the early universe thereby solving many of the problems of cosmology. Thus this scenario can successfully 
generate a homogeneous and isotropic Friedmann-like universe from initial conditions which, in the absence of inflation, 
would have resulted in a universe far removed from the one we live in today. However, still only a restricted set of 
initial data will lead to smooth enough conditions for the onset of inflation [120,203], so the issue of homogenization 
and isotropization is still not satisfactorily solved. Indeed, the initial conditions problem (that is, to explain why 
the universe is so isotropic and spatially homogeneous from generic initial conditions), is perhaps one of the central 
problems of modern theoretical cosmology. These issues were recently revisited in the context of brane cosmology 
[50,51]. 

A geometric formulation of the class of RS brane-world models was given in [208,212]. The dynamical equations 
on the 3-branc differ from the equations in GR. There are local (quadratic) energy-momentum corrections that are 
significant only at high energies and the dynamical equations reduce to the regular Einstein field equations of GR 
for late times. However, for very high energies (i.e, at early times), these additional energy momentum correction 
terms will play a critical role in the evolutionary dynamics of the brane-world models. In addition to the matter field 
corrections, there are nonlocal effects from the free gravitational field in the bulk, transmitted via the projection 
of the bulk Weyl tensor, that contribute further corrections to the Einstein equations on the brane. A particularly 
useful formulation of the governing equations has been derived [161,162] adopting a covariant approach [166,221]. 
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Recently, much effort has been devoted to understand the cosmological dynamics of RS brane-world mod- 
els. There are many generalizations of the original RS scenario which allow for matter with cosmologi- 
cal symmetry on the brane (Friedmann branes) [27,28,83] (in this case the bulk is Schwarzschild-Anti de 
Sitter space-time [30,184]), and non-empty bulks have also been considered including models with a dila- 
tonic scalar field in the bulk [10,168,175]. The Friedmann brane models have been extensively investigated 
[3,27,28,46,48,66,83,104,107,117,125,127,137,171,174,178,179,184,215,219] and inflationary scalar perturbations [165] 
and bulk perturbations [33,34,70,133,135,142,180,181] in these models have also been considered. The most important 
feature that distinguishes brane cosmology from the standard scenario is the fact that at high energies the Friedmann 
equation is modified by an extra term quadratic in the energy density. The dynamics of a brane-world universe filled 
with a perfect fluid have been intensively investigated [37,38,50,51]. It has been found there exist new regimes that 
are not inherent in the standard cosmology, such as stable oscillation [38] and the collapse of a flat universe [207]. 
Some features of brane-world inflation have been studied in [65,157-159,164,165] and the cosmological dynamics for 
scalar field models with an exponential potential have been described in [93,94,177]. 

The asymptotic dynamical evolution of spatially homogeneous brane-world cosmological models containing both a 
perfect fluid and a scalar field close to the initial singularity, where the energy density of the matter is larger than the 
brane tension and the behaviour deviates significantly from the classical GR case, has been studied in detail using 
dynamical systems techniques. The physical case of a perfect fluid with p = (7 — l)p satisfying 7 > 4/3, where 7 = 4/3 
corresponds to radiation and 7 = 2 corresponds to a massless scalar field close to the initial singularity, is of particular 
interest. We shall show that an isotropic singularity is a past-attractor in all orthogonal Bianchi models (including 
Bianchi type IX models), and is also a local past-attractor in a class of inhomogcneous brane-world models. The 
dynamics of these inhomogeneous brane-world models, together with the BKL conjectures [18-20], provide support 
for the conjecture that the initial cosmological singularity is isotropic in spatially inhomogeneous brane-world models. 
It is consequently argued that, unlike the situation in GR, it is plausible that the initial singularity is isotropic and 
that brane cosmology may naturally give rise to a set of initial data that provide the conditions for inflation to 
subsequently take place. 
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II. BRANE DYNAMICS 



The field equations induced on the brane lead to new terms that carry bulk effects onto the brane [208,212]: 

G/j,v — ^-9p.v "f" K + k S^ v £^ u , (1) 

where k 2 = 8tt/M 2 , k 2 = 8ir/Mp, and M p is the fundamental 5D Planck mass, which is typically much less than the 

effective Planck mass on the brane, the bulk cosmological constant A is negative and is the only 5D source of the 
gravitational field, and 

A = 6|J, A=±£ 2 (a+±k 2 A 2 ) . (2) 

We assume that A is chosen so that A = 0. The bulk corrections to the Einstein equations on the brane are of two 
forms: first, the matter fields contribute local quadratic energy-momentum corrections via the tensor S^, and second, 
there are nonlocal effects from the free gravitational field in the bulk, transmitted via the projection £^ v of the bulk 
Weyl tensor. The matter corrections are given by 

C _J_T^ ft T r T a J- — n rpaft (rp Q\2 / ry\ 

— 12 Q P- v 4 J ^ia J v T i\l)\xv •J- L a/3 1 y 1 a ) i \<J ) 

which reduces to 

S^v = ^p 2 u ll u l/ + ±p (p + 2p) , 

for a perfect fluid or minimally-coupled scalar field (or combination thereof) and is the 4-velocity comoving with 
matter and = g^ v + u^u v . The quadratic energy-momentum corrections to standard GR will be significant for 
K 4 p 2 > K 2 p (i.e., in the high-energy regime). 

Due to its symmetry properties, the projection of the bulk Weyl tensor can be irrcducibly decomposed with respect 
to u» [161,162] : 

—6 

K A 



where 



is an effective nonlocal energy density on the brane (which need not be positive), arising from the free gravitational 
field in the bulk. There is an effective nonlocal anisotropic stress 

on the brane, which carries gravitational wave effects of the free gravitational field in the bulk. The effective nonlocal 
energy flux on the brane is given by 

The local and nonlocal bulk modifications may be consolidated into an effective total energy-momentum tensor: 

G„„ = -Ag^ + k 2 T^ , (4) 

where 



— T^ v + „ £ M „ . (5) 



The effective total energy density, pressure, anisotropic stress and energy flux are [161,162] (for a comoving fluid 
the tilting case is treated later): 
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For an empty bulk, the brane energy-momentum tensor separately satisfies the conservation equations, V V T^ V = 0. 
The Bianchi identities on the brane then imply that the projected Weyl tensor obeys the constraint 



Consequently nonlocal bulk effects can be sourced by local bulk effects. The brane energy-momentum tensor and the 
consolidated effective energy-momentum tensor are both conserved separately. These conservation equations, as well 
as the brane field equations and Bianchi identities, are given in covariant form in [161,162]. 



In general, the 4 independent conservation equations determine 4 of the 9 independent components of £^ v on the 
brane. There is no evolution equation for the nonlocal anisotropic stress V^v Thus, in general, the projection of the 
5D field equations onto the brane does not lead to a closed system, since there are bulk degrees of freedom whose 
impact on the brane cannot be predicted by brane observers. These degrees of freedom could arise from propagating 
gravity waves in the bulk which are governed by off-brane 5D bulk dynamical equations. 

If the nonlocal anisotropic stress contribution from the bulk field vanishes, i.e., if V^ v = 0, then the evolution of 
£ M „ is fully determined. A special case of this arises when the brane is Robertson- Walker (RW) and in some special 
cases in which the source is perfect fluid matter and branes with isotropic 3-Ricci curvature (such as, for example, 
Bianchi I branes). Although these special cases can give consistent closure on the brane, there is no guarantee that 
the brane is embeddable in a regular bulk, unlike the case of a Friedmann brane whose symmetries imply (together 
with Z2 matching) that the bulk is Schwarzschild-AdS5 [30,184]. 

The nonlocal anisotropic stress terms enter into crucial dynamical equations, such as the Raychaudhuri equation 
and the shear propagation equation, and can lead to important changes from the GR case. The correction terms 
must be consistently derived from the higher-dimensional equations. In the analysis we shall primarily assume that 
the effective nonlocal anisotropic stress is zero (in the fluid comoving frame). This is supported by a dynamical 
analysis [57] and the fact that since corresponds to gravitational waves in higher-dimensions it is expected that 
the dynamics will not be affected significantly at early times close to the singularity [145,146]. This is the only 
assumption we shall make, and it is expected that inclusion of V^ v will not affect the qualitative dynamical features 
of the models close to the initial singularity. Indeed, recent analysis provides further support for V^ v = [55,56]. 
We expect that <~ Ug^ on dimensional grounds, and so for a Friedmann brane close to the initial singularity 
we expect that "P M „ ~ a 2 UC^ V (where is slowly varying), which is consistent with the linear (gravitational) 
perturbation analysis (in a pure AdS bulk background) [180]. Hence, V^ v is dynamically negligible close to the initial 
singularity. In addition, from an analysis of the evolution equation for Q M it can be shown that a small does not 
affect the dynamical evolution of U close to the initial singularity [50,102,161,162]. Let us discuss this in more detail. 
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A. Projected Weyl Tensor 
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1. Type N spacetimes 



In brane-world cosmology gravitational waves can propagate in the higher dimensions (i.e., in the 'bulk'). In some 
appropriate regimes the bulk gravitational waves may be approximated by plane waves [55] . For example, there might 
occur thermal radiation of bulk gravitons [141]. In particular, at sufficiently high energies particle interactions can 
produce 5D gravitons which are emitted into the bulk. Conversely, in models with a bulk black hole, there may be 
gravitational waves hitting the brane. At sufficiently large distances from the black hole these gravitational waves 
may be approximated as of type N [55] . Alternatively, if the brane has low energy initially, energy can be transferred 
onto the brane by bulk particles such as gravitions; an equilibrium is expected to set in once the brane energy density 
reaches a limiting value. We can study 5D gravitational waves that are algebraically special and of type N. 

Let us assume that the 5D bulk is algebraically special and of type N. This puts a constraint on the 5D Weyl tensor 
which makes it possible to deduce the form of the non-local stresses from a brane point of view. For a 5D type N 
spacetime there exists a frame £ a , £ a , m l a , i — 1, 2, 3 such that [58,59]: 

Cabcd = kCu\jt{ a rn\l c m 3 d y (11) 

The "electric part" of the Weyl tensor, £ a b = C ac bdn c n d , where n c is the normal vector on the brane, can for the type 
N bulk be written as 

Sab = C iaj [l a (mln c ) - ml(£ c n c )] [4(mjn c ) - mj(4n c )] , (12) 

where £ a bn b = £ a a = 0. Note that for a type N bulk we also have £ a b£ b = 0, which can be rewritten using the 
projection operator on the brane, as £ a b£ b = 0, £ b = g b c £ c - Hence, the vector ly, is the projection of the null vector 
£b onto the brane, and 

M b = -(IX) 2 . (13) 

In the case that £ a is time- like (£ a n a ^ 0) we can set parallel to £^. In this frame, the requirement £^£^, v = 
implies that we have U = = <2 M , and hence we can write £ llv = — {k/k^V^. This case was discussed in [55]. More 
importantly, in the case £ a n a = 0, in which £ fM is a null vector, £ ilv can be written 

£^ = -[^j eljv (14) 

Hence, this case is formally equivalent to the energy-momentum tensor of a null fluid or an extreme tilted perfect 
fluid. Using a covariant decomposition of £ t _ ll/ , the non-local energy terms are given by: 

U = e(£ v u v f 1 Q li = e{£ v u v )£ tl , V„ v = e£ { J v} . (15) 
The equations on the brane now close and the dynamical behaviour can be analysed [56]. Note that 

UVftv = QfiQu — 7^g fJ ,vQ\Q x , (16) 

so that in this determined completely by U and Q M . 

We can investigate the effect this type N bulk may have on the cosmological evolution of the brane. If we assume 
that we are in the regime of an isotropic past singularity and the cosmological evolution is dominated by an isotropic 
perfect fluid with equation of state p = (7 — l)p, the equations for U and <2 M are 



U + ^eU = Q, Q M + |eQ M = (17) 



(see also [141]). For the isotropic singularity, the expansion factor is given by = l/(-ft). Defining the expansion- 
normalised non-local density, U = U/& 2 and non-local energy flux, Q M = Q M /6 2 , we obtain 

17= 3^(37 -2)tf, Q M = 3^(37 -2)Q„. (18) 

Hence, the isotropic singularity is stable to the past with regards to these stresses if 7 > 2/3. Similarly, we have a 
Friedmann universe to the future with 9 = 2/(-yt), and thus 
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U= ^-(3 7 -4)C/, Qm = A(37-4)Q„. 



(19) 



This implies that if the isotropic fluid on the brane is less stiff than radiation at late times (7 < 4/3), then the 
Friedmann universe is stable to the future with respect to the non-local stresses. 

These qualitative results are further supported by a more detailed dynamical systems analysis [52,221] of the 
asymptotic behaviour of two tilting 7-law fluids in a class of Bianchi type VIo models [56] . In particular, the physically 
relevant case of interest here, namely that the second fluid is a null fluid or fluid with extreme tilt, was investigated. 
All of the equilibrium points were found and their stability determined, so that the local attractors were established. 
The effect a type N bulk may have on the cosmological evolution of the brane was then investigated. It was found 
that if 7 > 2/3, then the null fluid is not dynamically important at early times; that is, the effects of the projected 
Weyl tensor will not affect the dynamical behaviour close to the initial singularity. This supports the result that an 
isotropic singularity is a local stable past attractor. Also, if 7 < 4/3, the extreme tilting fluid is dynamically negligible 
to the future. Therefore, the null brane fluid is not dynamically important asymptotically at early times for all values 
of 7 of physical import, supporting the qualitative analysis above [56]. In particular, this implies that the effects of 
£ M „ are not dynamically important (at least in this class of models) in the asymptotic regime close to the singularity. 

2. Bulk spacetimes 

Most work on anisotropic brane-world universes use the effective 4D Einstein equations on the brane. Since the 
brane equations are not closed, various additional conditions on have been proposed in a rather ad hoc way. 
There has been some work on determining the effects of a non-zero on the behaviour of brane-world cosmological 
models [1,16,164,209,217]. For example, in [85] a bulk metric with a Kasner brane was presented. Since the Kasner 
metric is a solution of the 4D Einstein vacuum equations, the bulk metric is a simple warped extension; for the choice 
of bulk metric made in [85] it was shown that the brane can only be anisotropic if it contains a constant tension with 
p = P\ Other examples include the Schwarzschild black string solution [42,44], the spatially homogeneous Einstein 
space with a black hole [36] and the brane wave spacetimes [55], but all of these result in Bianchi models with an 
anisotropic fluid source. 

In order to have a fully consistent picture, we cannot avoid specifying the bulk geometry; i.e., we need to construct 
explicit anisotropic brane cosmologies by solving the full 5D vacuum Einstein equations with negative cosmological 
constant. This has non trivial implications for the brane itself because, in brane-world models, the Israel junction 
conditions (relating the extrinsic curvature to the matter on the brane) must be satisfied. For instance, if the brane 
contains perfect fluid matter, then the junction conditions will impose constraints on the components of the extrinsic 
curvature. 

In [39] a moving brane admitting a spatially homogeneous but anisotropic Bianchi I 3D slicing in a static anisotropic 
bulk was considered. Having solved the bulk Einstein equations, it was shown that the junction conditions induce 
anisotropic stresses in the matter on the brane (i.e., the brane cannot contain a perfect fluid). Indeed, from the bulk 
Einstein equations and the junction conditions, it follows that this anisotropic stress can only vanish if the bulk is 
isotropic and hence the brane is isotropic. It was also noted that the anisotropic stresses obtained on the brane in 
[39] are fixed by the bulk metric. Some explicit Bianchi I brane-world cosmological models using the freedom still 
available in embedding the 3-brane were constructed. In one example the cosmological models do not isotropize as the 
initial singularity is approached, but in this example w — ► —1 as t — > 0! (In other examples, the models do isotropize 
in the past, but the matter content never behaves as a perfect fluid.) 

This work was generalized in [82] to a class of non-static bulks. The bulk was assumed to be empty but endowed 
with a negative cosmological constant and explicit analytic bulk solutions with anisotropic 3D spatial slices were 
found. The Z 2 -symmetric branes were embedded in the anisotropic spacetimes and the constraints on the brane 
energy-momentum tensor due to the 5D anisotropic geometry were discussed. The question of whether it is possible 
to find a bulk geometry and a brane with perfect fluid matter consistent with the Israel junction conditions was 
studied. Einstein's equations were integrated explicitly in the cases in which the shear does not depend on the extra- 
dimension and in which the metric is separable. Remarkably, in the separable bulk solution the analysis implies that 
the brane can support perfect fluid type matter for certain parameter values. In a particular example, by solving for 
the cosmological evolution on the brane it was found that as r — > —00, p converges towards the RS tension and the 
effective pressure is negative (but converges towards zero more rapidly than the effective energy). Therefore, unlike 
[39], it was found [82] that it is possible to find a perfect fluid anisotropic brane in some such bulk geometries, but 
again there is no flexibility in the choice of the perfect fluid in contrast with the isotropic case (i.e., only a very 
particular type of perfect fluid is compatible with the given bulk geometry) . 
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Therefore, it does seem difficult to construct examples of anisotropic branes with perfect fluid matter. However, 
these solutions [39,82] are very special examples and say nothing about the general case. An exact solution in the 
bulk gives rise to a precise relationship between the density and the pressure on the brane. But exact closed form 
solutions are few sparse; indeed, even in 4D GR there are only a few very special exact Bianchi solutions (models are 
usually only specified up to a set of ODE), so this approach will not shed light on the general case. Moreover, the 
isotropic singularity conjecture is not contradicted by these results. First, it is not necessary for V^v to be zero but 
rather that it be dynamically negligible (i.e., some appropriate normalized quantity is negligible) as the singularity 
is approached. Second, the values for the parameter 7 is in the wrong range in the examples of [39] and [82] (the 
isotropic singularity theorems only apply for 7 > 4/3). 

In order to illustrate this point further, let us consider the following Bianchi I metric: 

ds 2 = -£2(7-1) [l + a f7-2]-T^ dt 2 +t 4/ 3[1 + af f-2 ] - 2< \ + { ;-; ) ^2 

_ 2(l-co.i»-V3.iml) _ 2(l-cosi/. + y / 3aini/.) 

+t 4 / 3 [l + af- 2 ] 5*7=35 dy 2 + t 4 / 3 [l + crf^ 2 ] dz 2 (20) 

where a and tp are constants and 7 > 2. For a = 0, the metric is the Binetruy, Deffayet and Langlois or BRW 
(see later) solution with a(T) = [27,28,83] (where T e (t). We note that for the metric (20) G/ = diag 
(-A,B,B,B), where B = (7- I) A and 

A = A r 2 ~[l + at~~ 2 ]W^. 

We shall show that a (anisotropic) Bianchi I perfect fluid RS brane can be embedded in a 5D Einstein space using 
the Campbell-Magaard theorem. Following Dahlia and Romero [69] we consider the 5D metric in Gaussian normal 
form (with 5th coordinate tj. If we wish to embed a RS brane into a Z 2 symmetric bulk (with cosmologial constant 
A) we need to consider the equations (on the brane t = £ , a constant). 

T^ = p+(p + p)Q = 0, (21) 



-{T^ - -T 2 ) = — [p 2 + 3pp + p\- 3Ap - 2A 2 ] = 2A + R, (22) 

where is the energy- momentum tensor for a perfect fluid and R is the scalar curvature of the hypersurface I = £q. 
If a solution of Eqs. (21) and (22) can be found, then an analytic solution of the 5D bulk equations = —Ag^ 
(satisfying the Eqs. (21) and (22) on the brane) is guaranteed by the Cauchy-Kowalewski theorem [69]. 
By a direct calculation, for metric (20) we obtain 

R = R t- 2 ~[l + at~- 2 }^ (23) 
6 = 9 t-~[l + at~- 2 ]^[l + ^at~- 2 }. (24) 

With these expressions, Eqs. (21) and (22) can always be solved to find p = p(t),p = p(t). To lowest order in 
t, we have that £0 = i?o^ 2 ~(f + 0(t)), 6 = 9 t^~(l + G(t)), and we can solve Eqs. (21) and (22) to obtain 
/i = p = pot— 7(1 + 0{t)) and p — (7 — l)p(l + 0(t)). That is, to lowest order we obtain the BRW solution. 

Therefore, using extended Campbell-Magaard theory [69], we have demonstrated that there are exact anisotropic 
(of Bianchi type I) perfect fluid brane cosmological models that can be embedded into a 5D bulk Einstein space (see 
also [128]). This construction does not, of course, gaurentee that the bulk is regular. Note that for the example (20) 
constructed, as t — > the model isotropizes (towards the BRW solution) and the singularity is isotropic. 



3. Dynamical systems approach 

In order to illustrate the possible effects of on the dynamics of a brane-world cosmological model, an alternative 
approach was taken in [110]. Again it was assumed that the effective cosmological constant on a Bianchi I brane is 
zero, and the matter on the brane is of the form of a non-tilted perfect fluid. Utilizing new dependent dynamical 
variables, the effective Einstein field equations yield a dynamical system describing the evolution of the Bianchi I 
brane-world. In this analysis there is no evolution equation for the quantity 
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~ k 2 X aH 2 



(25) 



modelling the non-local bulk corrections. In order to obtain a closed system of equations, it was then assumed that V 
is a diffcrcntiable function of the remaining dynamical variables; that is, V = V(T,, 0, fl\, flu). With this assumption, 
the resulting brane- world spacetimes represented by the equilibrium points of the system are geometrically self-similar 
[40]. 

A number of models were investigated. In general, the dynamics in the case V = Vo (a constant) is quite different 
from that found currently in the literature. For P^l (Vo ^ 0), none of the equilibrium points are isotropic, and all 
have some form of dark matter with U ^ 0. It appears that the non-local bulk corrections have a significant impact on 
the dynamical behaviour and the initial singularity is not isotropic. In the case V = "PsE" (a > 0), the BRW brane- 
world solution Tb may not represent the past asymptotic state (there are parameter values of Ps such that the initial 
singularity is represented by a Kasner-like solution). In the case V = Vn u n u (J3 > 0) [16,209], T b is always a source, 
although there will exist cosmological models that will asymptote to the past towards anisotropic bulk-dominated 
Kasner-like models (i.e., the initial singularity is not necessarily isotropic). Thus in each of these three cases, the 
effects of the non-local "P M „ could be significant. Unlike the standard cosmological brane- world situation in which Tb 
represents the 'generic' initial behaviour in brane-world cosmologies containing ordinary matter, it was found that 
Kasner-like bulk-dominated models can also act as sources and that Tb does not even exist as an equilibrium state in 
one of the cases studied. The future behaviour of the models is determined primarily by the parameter 7, and was 
studied in [110]. 

We note that, from (25), V is not defined in the (isotropic limit) as a — > 0. This is problematic when attempting to 
analyse isotropic equilibrium points; indeed, the isotropic limit cannot be studied and the analysis in [110] is essentially 
concerned with the local stability of non-isotropic equilibrium points. However, the DEs may be well-defined and have 
an isotropic limit depending on the chosen ansatz for V . In the second and third cases, a well-defined isotropic limit 
exists and in these cases Tb is a local attractor. In the first case (V constant), there is no isotropic limit possible. 
Indeed, the models do not even allow for an isotropic solution at any time (the BRW models do not exist), and 
any conclusions concerning an isotopic singularity are inappropriate in this case. Consequently the analysis in [110] 
supports the existence of a local isotropic singularity. 
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III. DYNAMICAL EQUATIONS 



The general form of the brane energy-momentum tensor for any matter fields of perfect fluid form (including scalar 
fields) can be covariantly decomposed. We adopt the notation of [161,162]. Angled brackets denote the projected, 
symmetric and tracefree part: 

V (li) = h„ v V v , W {IU>) = [h^Kf - §fc°"V] W «P . (26) 

with round brackets denoting symmetrization. An overdot denotes u u V v (i.e., timelike directional derivative), O = 
V^u^ is the volume expansion rate of the congruence, = = Ai^ is its 4-acceleration, <j^ v = D/ ll u v \ is its 
shear rate, and ui^ — — icurlw^ = is its vorticity rate. The covariant spatial curl is given by [161,162] 

curl^= W D a y^, cmlW^=s a ^B a W^ v) , (27) 

where e^ va is the projection orthogonal to of the brane alternating tensor, and is the projected part of the 
brane covariant derivative (i.e., spatial derivative), defined by 

Djuf" '•••/3 = (V M F a -... /J ) ± = h l /h a 1 ■ ■ ■ h p s V v FV"... s . (28) 

The remaining covariant equations on the brane are the propagation and constraint equations for the kinematic 
quantities <d, A^, lo^ 1 a^ v , and for the nonlocal gravitational field on the brane. The nonlocal gravitational field on 
the brane is given by the brane Weyl tensor C^ ua fj. This splits into the gravito-electric and gravito-magnetic fields 
on the brane: 

— C t _ iai/ f}u a u t3 = 2?( M „) , = \e ^ a f}C al3 ^u 1 = . (29) 

The Ricci identity for and the Bianchi identities V^C^a/j — V[ jU (— R v ] a + \Rgv]a) give rise to the evolution and 
constraint equations governing the above covariant quantities once the Ricci tensor i? M „ is replaced by the effective 
total energy-momentum tensor from the field equations [75,77,163]. In the following Rh v is the Ricci tensor for 
3-surfaces orthogonal to on the brane and R^- = h^ v R^ v . Note that in the GR limit, A -1 — > 0, we have that 
£„„-0. 

4- Vorticity-free models with a linear barotropic equation of state 

In the case of a perfect fluid with an equation of state p = (7 — l)p, we obtain the following equations when the 
vorticity is zero (A ^ 0) [161,162]: 

A, = -^1-B^p, (30) 
IP 

P = -7&P, (31) 

U + §OW + D tl Q^ + 2A» Q„ + a^Vpv = , (32) 

Q(n) + feQ^ + §D„W + \UA^ + D"P„ V + A V V^ + a^Q" = -^jpB^p. (33) 



6 + ±9 2 + a^ v - DM„ + A^ + i« 2 (3 7 — 2)p — A = -f^(3 7 - l)p 2 - ~^rU . (34) 

K A 

3 

V(nv) + 3©<V, + E^ v - T)^A V ) + a a ^a v ) a - A^A V ) = —^-V^ ■ (35) 



K 2 X 



E(iiv) + - curltf^ + ^po> - 2A a e aP ^H u) f3 - 2,a a{lx E v 



2A P u » v K 2\ 



4Ua^ + 37V) + ©?V + 3D {ll Q v) + 6A {li Q v) + Za a {l J> v)a . (36) 



H {ltv) + @H„ V + cmlE„ v - ZcT a ^H v) a + 2A a e a ^E v / = 4y [ curl ?V + ^ a {^ v)a pQ t3 } . (37) 

K A 
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D^-fD^-J^. (38) 
curl a^- if ^ =0. (39) 



D"^ - i« 2 D M p - [a, H}„ = -—^-D^p + — (2D M W - 26Q M - 3D"P„ V + 2>a lxv Q v ), (40) 

oA /v A 



'- 2 D^-[a,fl] M = ^D^-f ' 

1 



K 2 A 

and the Gauss-Codazzi equations on the brane: 



D"ff M „ + [a, E]„ = -— (3curl Q M + 3[<x, V]„) , (41) 



^ + |e 2 - - 2 K 2 p - 2A = ^p 2 + . (43) 



5. Integrability conditions for V^ v — 0, Q M = 
When "P M „ = and Q M = 0, these equations simplify. Using Eq.(30) , Eq. (33) becomes for 7 7^ 1: 



4(7 -l)U K 4 - 
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(44) 



Taking the directional time derivative of Eq. (44), and using the relations for interchanging space and time derivatives, 
e-g-, 



(45) 



we obtain the integrability condition 

"4 



-(3 7 -4)W-^V 



4 

K 7 



— (3 7 - 4) (7 - \)U + ^(37- 1)P 2 



Taking the directional time derivative of (38) and interchanging space and time derivatives, we obtain 

^I^[12W + kV]-D„p = 0. 
7 p 



(46) 



(47) 



Hence for 7 7^ 1 and [12U + K 4 p 2 } 7^ 0, we have 



so that 



D M p = 0, 



= D^U = D M 9 = D"V„„ = D m( j 2 = 

and hence, in general, the brane is spatially homogeneous. The special case 7 = |, U + K 4 p 2 = 0, in which the 
integrability conditions yield no constraints, corresponds to the case in which there are no corrections to the general 
relativistic equations. Finally, in the special case 7 = 1, in which p = and = 0, we obtain analogues of Eqs (44) 
and (46) from which it follows that p = 0, and hence 6 = 0. 
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A. Spatially Homogeneous Branes 

In the case of spatial homogeneity we have that 

w„ - 0, Aft = 0, D^G = D M p = D^p = B^U = 0. 
With a general equation of state, we then obtain 

p + Q( P + p) = 0, 



(48) 



U + | BU + D"Q M + a^V^ = , 

e< M > + feQ^ + D "?V + = o. 



(49) 
(50) 



6 + ie 2 + <7 M „ + \n 2 {p + 3p) - A = -£(2p + 3p)p - ^-U . 

K A 



(51) 
(52) 



+ QE^ - curltf^ + i K 2 (p +p)a liV - 2>u a{ll E v) a 

A K A L 



# <At „) + + cuilE^ - "ia a{ll H u) a = — [curl TV + a a ( M £ I ,) Q(3 S /3 ] 

K. A 



(53) 
(54) 
(55) 



curla^y - = . 



D"^ - [a, = - — (29Q M + 3D"7V - 3^<2"), 

A 



D"^ + [a, £] M = — (3curl Q M + 3[a, V}„) . 

K A 



The Gauss-Codazzi equations on the brane are: 



R 1 - + |6 2 - - 2n 2 p - 2A = ^p 2 + ^Ui . 

A hi A 



(56) 
(57) 

(58) 

(59) 
(60) 



1. RW brane 



It follows from the symmetries that for a RW brane Q M = = V^. The generalized Friedmann equation (60) on a 
spatially homogeneous and isotropic brane is [27,28,83] 



^V(i + £Ha-* + §(^)\ 



(61) 



12 



where k = 0, ±1. The nonlocal term that arises from bulk Coulomb effects (also called the dark radiation term) is 
strongly limited by conventional nucleosynthesis [145,146]: (U / p) nuc \/ k 2 X < 0.005, where U = U (a /a) 4: . A more 
stringent constraint comes from small scale gravity experiments [141]. This implies that the unconventional evolution 
in the high energy regime ends, resulting in a transition to the conventional regime, and the nonlocal term is sub- 
dominant during the radiation era and rapidly becomes negligible thereafter. The generalized Raychaudhuri equation 
(51) (with A = = U) reduces to 



H + H 2 = -\n 2 



p + 3p+ (2p + 3p)^l . (62) 
A. 



The bulk metric for a flat Friedmann brane is given explicitly in [27,28,83]. In natural static and spherically 
symmetric coordinates, the bulk metric is Schwarzchild-anti de Sitter (with the .^-symmetric RW brane at the 
boundary) and is given in [30,184]. The RW brane moves radially along the 5th dimension, with R — a(T), where a 
is the RW scale factor, and the junction conditions determine the velocity via the Friedmann equation for a. Thus 
one can interpret the expansion of the universe as motion of the brane through the static bulk. In Gaussian normal 
coordinates, the brane is fixed but the bulk metric is not manifestly static. The bulk black hole gives rise to dark 
radiation on the brane via its Coulomb effect. The mass parameter of the black hole in the bulk is proportional to 
U . In order to avoid a naked singularity, we assume that the black hole mass is non- negative. 

In general, U ^ in the Friedmann background [27,28]. Thus for a perturbed Friedmann model, the nonlocal bulk 
effects are covariantly and gauge- invariantly described by the first-order quantities D M W, <2 M , V^. 



B. Static Branes 



Assuming spatial homogeneity (so that = 0) with Q M = Vy, v = 0, in the isotropic curvature subcase with 
Rf^v) — 0, so that R^- = ka~ 2 and 

J?- 1 = -fei?- 1 , (63) 

we have that 



P = -1&P, (64) 
U = -\m, (65) 

6 + ie 2 + 2<7 2 + i K 2 (3 7 - 2)p - A = -fx(37 - i)p 2 - 4t w • ( 66 ) 

K A 

i? ± + |e 2 -2a 2 -2« 2 p-2A = ^p 2 + ^-W, (67) 

A k A 

(j 2 = -26cr 2 , (68) 

and the system of equations is closed. 

In the static models we have that 6 = (and we shall assume that a = 0), and 

± K 2 (3 7 - 2) Po - A = - £(37 - 1)pI - 4t"o • (69) 

K A 

< + -2^ - 2A = ^-p 2 + -^-U a , (70) 

A K A 



where po, Uq and R ± o are constants. These constitute two algebraic equations for four constants (and 7, which need 
not be equal to 2/3). Since Uq need not be positive, i? ± o need not be positive definite. There are static RW branes 
that can be embedded in a bulk that is not Schwarzschild-AdS [91]. 



C. The EGS Theorem in Brane- World Models 



The isotropy of the cosmic microwave background (CMB) radiation has crucial implications for the spatial homo- 
geneity of the universe. If all fundamental observers after last scattering observe an isotropic CMB, then it follows 
from a theorem of Ehlers, Geren and Sachs (hereafter EGS) [74] that in GR the universe must have a RW geometry. 
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More precisely, the EGS theorem states that if all observers in a dust universe see an 'isotropic radiation field', which 
is implicitly identified with the CMB, then that spacetime is spatially homogeneous and isotropic. This can trivially 
be generalised to the case of a geodesic and barotropic perfect fluid [78,79]. However, as has been emphasised recently 
[13,47], the resulting spacetime will be RW only if the matter content is that of a perfect fluid form and the observers 
are geodesic and irrotational. The EGS theorem has recently been investigated in inhomogeneous universes models 
with no n- geodesic observers, and inhomogeneous spacetimes have been found which also allow every observer to see 
an isotropic CMB [13,47]. 

The EGS theorem is based on the collisionlcss Boltzmann equation and on the dynamical field equations. Bulk 
effects in brane- world models do not change the Boltzmann equation, but they do change the dynamical field equations 
[49,161,184,208,212], so that in principle the resulting geometry need not be RW on the brane. The theorem is 
concerned with universes in which all observers see an isotropic radiation field. It follows from the multipole expansions 
of the Einstein-Boltzman equations for photons in a curved spacetime that even in the brane- world scenario a spacetime 
with an isotropic radiation field must have a velocity field of the photons which is shearfree and obeys [78,79] 

= D M g, 6 = 34, (71) 

where q is a function of the energy density of the radiation field. Any observers traveling on this congruence will then 
observe the isotropic radiation. This velocity field is also a conformal Killing vector of the spacetime, and hence a 
spacetime admitting an isotropic radiation field must be conformally stationary, in which case the metric can be given 
in local coordinates by 

ds 2 = e 2 9(*.* a ) {-dt 2 + h aP dx a dx (3 } (72) 

where h a p{x' 1 ) can be diagonalised. If q = q(t), then the acceleration is zero. 

Therefore, assuming (71) and a = (as well as w = 0), we immediately have that 



v„e = 4x Q ^ (73) 



H^ = 0. (74) 
From Eqs. (30) and (31) we then obtain (assuming p = (7 — \)p) 

D M e - m^q) = = 3(7 - i)D„e , (75) 

which implies that (for 7 ^ |) 

07,6 = 0. (76) 

(When 7 = I we obtain p = po<Z~ 4 )- Consequently, from Eq. (73) we have that 

Qij, = . (77) 

From (75) we have that D M (4) = and [D v q\ = 0, which then implies that = 0. This can best be seen by using 
local coordinates (72), in which these conditions yield q{t,x a ) = f(t) + gi(x a ) and e~ q = e~f^ + g2(x a ), so that for 
consistency either f(t) = (whence 9 = 0) or q — q(t), and hence A^ = 0. Therefore, we find that 

A„ = 0, (78) 

and hence D^p = 0. 

Taking = 0, = 0, D^p = 0, D^B = 0, A^ = 0, which implies that 



3 



Eqs. (33) and (40) then yield 

D M W = 0, D"7V = 0, D V £7 M „ = , (80) 

and Eqs. (36) and (42) become 
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2 ± 

E (^u) + t^iE^ = , R^ v) = lEpv . (81) 
For the metric (72) with q = q(t), it then follows that 

RU=RU(x T ), (82) 
where R^ v is the Ricci tensor of the 3-metric h a p{x 1 ) [60]. Consequently 

dt[Rt v) }=0, (83) 

and hence dt[E [_,,„] = and dt[P^v\ = 0. The remaining non-trivial equations are the conservation law for U, the 
generalized Raychaudhuri equation, and the generalized Friedmann equation. 

1. Discussion 

The 5D metric in the bulk is given in Gaussian normal coordinates by 

ds 2 = g AB dx A dx B = g ab (t, x'< ,y)dx a dx b + dy 2 , (84) 

where A, B = 0, 1 — 3, 4 = y and a, b = 0, 1 — 3, the brane is located at y = and the unit normal to the brane 
is n A — Sy. Assuming that the metric functions g a b are smooth, they can be expanded in powers of y close to the 
brane up to 0(y 2 ). The metric functions are determined up to O(y ) by Eq. (72). The extrinsic curvature on the 
brane, K\ B , which is discontinuous at y = 0, restricts the form of g a b(t, x 1 , y) to 0(y). Assuming that there arc no 
5D fluxes the stress-energy of the bulk, Tab, in the neighbourhood of the brane satisfies T 0a = 0, which is valid for 
scalar fields, perfect fluids, a cosmological constant and combinations thereof, and we have that go a — (which then 
implies that <2 M = 0). The non-trivial components of g a \, up to 0(y 2 ) can then be written (with a slight change of 
notation) in the y > ( + -region) as 

g 00 = -(l + f(t)y + F(t,x->)y 2 ), (85) 

g a p = Q(t)(l + q(t)y)h af} + G af3 (t, x->)y 2 . (86) 
Noting that £ a b = C a cbDn c n D — C a 5b5, by a direct calculation we find that on the brane (y = 0) 

do5 ° 5 = wtf [t) + wkw) {QQ{t)F{t ' xl) ~ R[xl) - 2G ^ X ^ ( 8? ) 

C a505 = -^R a(j {x-<) + ^( 3 R(x->)+T(t) + 2Q(t)F(t,x-*) + 2G(t,x-y))h a p(x-y) - ^G a p(t,x^), (88) 



where 



T(t) = Q tt - l(g t ) 2 - \Qfq + Qq 2 - \Q!\ (89) 
G{t,xi) = h a ^G aP , (90) 



and C M 505 = (so that Q M = 0). 

From £{ f _ LV ) = —Rf^ix 1 ) an d £oo — — (t) , we then find that 



and 



G aP (t,x<) = RifAx 1 ) + \(G(t,x>) - R^{x^))h a ^), (91) 



G(t, x<) = |.F(t) + + 3 ^(W> x ^ - \ R ^^ T )- ( 92 ) 
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Finally, for a bulk source consisting of a combination of scalar fields, perfect fluids and a cosmological constant, we 
then also have that 

F(t, x r ) = F(t), G(t, x' 1 ) = G(t), (93) 

and it follows that i?- L (x 7 ) is constant. 

We can solve the full equations in the bulk iteratively. In general it seems that the remaining bulk field equations 
consist of five ODEs for the five non-trivial functions (of t) in the metric and additional matter fields and so there 
will be many possible solutions. It may be possible to integrate the bulk equations explicitly in the case of a 5D 
cosmological constant. 
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IV. ESSENTIAL FEATURES 



The generalized Friedmann equation, which determines the volume expansion of the universe or the Hubble function 
H = |0, in the case of spatially homogeneous cosmological models is 

a, -5* , '( 1 + s)-5« l + 5°' + 5 A + S- <94 » 

where 2a 2 = <J ab <J a b is the shear scalar and H = a/ a [161,208,212], and the background induced metric satisfies 

D^ = Q^ = 7V = 0. (95) 
For a flat isotropic brane (with A = = U) we obtain 

The conservation equation is given by 

p + 3H(p + p) = 0. (97) 
For a minimally coupled scalar field the energy density and pressure are, respectively, 

P=\tf + V{4>), P =^-V(cb), (98) 

so that the conservation law is equivalent to the Klein-Gordon equation 

4> + 3H<p + V , (<p) = 0. (99) 

At early times, assuming that the energy density is dominated by a scalar field and that p is very large (and U = 0), 
the conservation equation (97) becomes 

p = -3</> 2 p, (100) 

so that p is monotonically decreasing and the models will eventually evolve to the low density regime and hence the 
usual general relativitistic scenario. 



A. Anisotropic and Curved Branes 



There are many reasons to consider cosmological models that are more general than RW, both spatially homogeneous 
and anisotropic and spatially inhomogeneous. The 3-curvature in RW models is given by R 1 ^ = % where k = 0,±1 
is the curvature constant. The structure of the initial singularity in RW brane- models was studied in [119]. An 
equivalent 3-curvature occurs in orthogonal spatially homogeneous and isotropic curvature models, and a similar term 
occurs in other models such as Bianchi V cosmological models. For a RW model on the brane, Eq. (95) follows, which 
implies that U = U(t). In the case of Bianchi type I we get Q p = but we do not get any restriction on 7"** '. Since 
there is no way of fixing the dynamics of this tensor we can study the particular case in which it is zero. Thus, in RW 
and Bianchi I models the evolution equation for U is U = —4HU [161,162], which integrates to U = ^r, which has 
the structure of a radiation fluid (sometimes referred to as the "dark radiation term", but where Uq can be negative). 

A Bianchi I brane is covariantly characterized in [164]. The conservation equations reduce to the conservation 
equation (97), an evolution equation for the effective nonlocal energy density on the brane U, and a differential 
constraint on the effective nonlocal anisotropic stress V^v In the Bianchi I case, the presence of the nonlocal bulk 
tensor V^v m the Gauss-Codazzi equations on the brane means that we cannot simply integrate to find the shear as in 
GR. However, when the nonlocal energy density vanishes or is negligible, i.e., U = 0, then the conservation equations 
imply o^V^v = 0. This consistency condition implies a^V^ = QV^V^vj n 2 \. Since there is no evolution equation 
for V^v on the brane [161,162], this is consistent on the brane. This assumption is often made in the case of RW 
branes, and in that case it leads to a conformally flat bulk geometry [27,28,83]. The shear evolution equation may be 
integrated after contracting it with the shear, to give 
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^^ = -6^ (101) 



where So is constant. Bianchi I models on the brane have been studied with a massive scalar field [164], with a 
perfect fluid with linear equation of state [37,38,217] and a perfect fluid and scalar field [186]. A similar shear term 
occurs in Bianchi cosmological models (such as Bianchi type V models) in which the hypersurfaces of homogeneity 
are orthogonal to the fluid velocity. 

The conservation equation is given by equation (97). If we assume that the matter content is equivalent to that of 
a non-tilting perfect fluid with a linear barotropic equation of state (i.e., p = (7 — l)p, where p > 0, and 7 £ [0, 2]) 
the conservation equation then yields p = po a ~ 37 , where po > 0. A dynamical analysis of scalar field models indicates 
that at early times the scalar field is effectively massless. A massless scalar field is equivalent to such a perfect fluid 
with a stiff equation of state parameter 7 = 2. 

We can therefore write down a phcnomenological (spatially homogeneous) generalized Friedmann equation: 

H -3^ + 6X^-^ + 3 A+ ^ + ^- (102) 

In many applications the 4D cosmological constant is assumed to be zero [193,194]; here we shall assume that if it is 
non-zero it is positive; i.e., A > 0. 

In particular, the cosmological evolution of the RS brane-world scenarios in RW and the Bianchi I and V perfect 
fluid models with a linear barotropic equation of state with U = were studied in [37,38] using dynamical systems 
techniques. This work was generalized in [38], in which the 5D Weyl tensor has a non- vanishing projection onto 
the three-brane where matter fields are confined, i.e., U 7^ 0; in the case of RW models the study was completely 
general whereas in the Bianchi type I case the Weyl tensor components were neglected (the Bianchi V models were 
not considered). 

The particular models considered so far, in which the curvature and shear are given by the expressions in the above 
phcnomenological equation, are very special. In particular, the Bianchi I and V models are not generic, and so the 
study of the dynamics of these models does not shed light on the typical behaviour of spatially homogeneous brane 
models. Clearly, more general scenarios must be considered to obtain physical insights, and we shall study Bianchi 
type IX models in detail later. 



1. Late times: 



Recent observations of distant supernovae and galaxy clusters seem to suggest that our universe is presently un- 
dergoing a phase of accelerated expansion [88,190,199,210], indicating the dominance of dark energy with negative 
pressure in our present universe. The idea that a slowly rolling scalar field provides the dominant contribution to the 
present energy density has gained prominence in recent times [195,220,223]. In [173] a potential consisting of a general 
combination of two exponential terms was considered within the context of brane cosmology to determine whether it is 
possible to obtain both early time inflation and accelerated expansion during the present epoch through the dynamics 
of the same scalar field. Such a potential, motivated by phenomenological considerations from scalar tensor theories 
of gravity and higher dimensional quantum effects or a cosmological constant together with a 4D potential, has been 
recently claimed to conform to all the current observational constraints on quintessence (for certain values of the 
parameters) [12]. The brane-world dynamics was assumed to be governed by the modified Friedmann Eq. (96) where, 
for reasons of simplicity, the contributions from bulk gravitons and a higher dimensional cosmological constant were 
set to zero [65,116,165]. It was found that the brane-world inflationary scenario is feasible with steep potentials (in 
contrast to the situation in standard cosmology). During inflation the necessity of sufficient inflation and the COBE 
normalised amplitude of density perturbations can be used to fix the values of the brane tension and the scale of the 
potential at this stage. The emergence of a second phase of accelerated expansion necessitates the introduction of a 
second exponential term in the potential. A successful brane quintessence model with a scalar field (and a power-law 
potential) and a dark radiation component (Uo ^ 0) was discussed in [10,168,175], with particular emphasis on how 
the quadratic matter term affects the evolution of the scalar field. 

When a positive cosmological constant is present (A > 0), the de Sitter model is always the global attractor for 
U > 0. For U < 0, models can (re)collapse (even without a positive curvature); i.e., there exist (rc)collapsing models 
both for RW and Bianchi type I models for any values of 7. This means that for U < the de Sitter model is 
only a local attractor [37,38] (and the cosmic no- hair theorem is consequently violated). In [207] a set of sufficient 
conditions which must be satisfied by the brane matter and bulk metric so that a homogeneous and anisotropic brane 
asymptotically evolves to a de Sitter spacetime in the presence of a positive cosmological constant on the brane was 
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derived. It was shown that from violations of these sufficient conditions a negative nonlocal energy density or the 
presence of strong anisotropic stress (i.e., a magnetic field) may lead the brane to collapse. For U < and positive 
curvature (as in the k = 1 RW models), not only can models recollapse, but there exist oscillating universes in which 
the physical variables oscillate periodically between a minimum and a maximum value without reaching any spacclikc 
singularity [38]. Indeed, when there are bulk effects, such as for example in case of a negative contribution of the 
dark-energy density term U, the singularity theorems are violated and a singularity can be avoided. 

2. Intermediate times: inflation 

The qualitative properties of Bianchi I brane models with a perfect fluid and a linear equation of state (U = and 
hi 7^ 0) were studied in [37,38]. A variety of intermediate behaviours can occur, and the various bifurcations in these 
models were discussed in some detail. In addition, the maximum value of the shear in Bianchi I models was studied 
in [217]. In particular, models with a positive curvature can recollapse, although the condition for recollapse in the 
brane- world can be different to that in GR [37,38]. 

Recent measurements of the power spectrum of the CMB anisotropy [9,25,140] provide strong justification for 
inflation [150,160,153]. Since the Friedmann equation is modified by an extra term quadratic in energy density p at 
high energies, it is of interest to study the implications of such a modification on the inflationary paradigm. The issue 
of inflation on the brane was first investigated in [165], where it was shown that on a RW brane in 5D anti de Sitter 
space, extra-dimensional effects are conducive to the advent of inflation. It has also been realized that the brane-world 
scenario is more suitable for inflation with steep potentials because the quadratic term in p increases friction in the 
inflaton field equation [165]. This feature has been exploited to construct inflationary models using both large inverse 
power law [116] as well as steep exponential [65] potentials for the scalar field. 

In particular, for a RW brane inflation at high energies (p > A) proceeds at a higher rate than the corresponding 
rate in GR. As noted earlier, this introduces important changes to the dynamics of the early universe, and accounts 
for an increase in the amplitude of scalar [165] and tensor [145,146] fluctuations at Hubblc-crossing, and for a change 
to the evolution of large-scale density perturbations during inflation [102]. The condition that radiation domination 
sets in before nucleosynthesis can be used to impose constraints on the parameters of brane inflation models [65,116]. 
A definite prediction of these models is the parameter independence of the spectral index of scalar density pertur- 
bations [65,116,165]. The behaviour of an anisotropic brane-world in the presence of inflationary scalar fields was 
examined in [164]. The evolution equations on the brane (with So ^ 0) for a minimally coupled massive scalar field 
with V = ^m 2 (f) 2 was studied and it was shown that, contrary to expectations, a large anisotropy does not adversely 
affect inflation [164]. In fact, a large initial anisotropy introduces more damping into the scalar field equation of 
motion, resulting in greater inflation. After the kinetic term and the anisotropy have dropped to a very small value, 
a transient anisotropy dominated regime begins and the anisotropy then rises, gradually approaching its asymptotic 
'slow-roll' value for the simple 'chaotic' potential. The kinetic energy does not remain constant but gradually increases 
as the field amplitude decreases during slow-roll. Since the decay of anisotropy is found to be gcncrically accompa- 
nied by a corresponding decrease in the kinetic energy of the scalar field [164], this effect leads to greater inflation. 
Subsequently, the anisotropy dissipates. 

In RS brane-worlds, where the bulk has only a vacuum energy, inflation on the brane must be driven by a 4D scalar 
field trapped on the brane [152]. In more general brane-worlds, where the bulk contains a 5D scalar field, it is possible 
that the 5D field induces inflation on the brane via its effective projection [8,32,45,86,107,108,132,136,143,170,176]. 
More exotic possibilities arise from the interaction between two branes, including possible collisions, which is me- 
diated by a 5D scalar field and which can induce either inflation [71,126] or a hot big-bang radiation era, as in 
the "ekpyrotic" or cyclic scenario [124,129,214], or colliding bubble scenarios [29,89,90] (for colliding branes in an 
M theory approach see [17,64]). In general, high-energy brane-world modifications to the dynamics of inflation on 
the brane have been investigated [21,22,105,172,174,177,200]. Essentially, the high-energy corrections provide in- 
creased Hubble damping making slow-roll inflation possible even for potentials that would be too steep in standard 
cosmology [61,65,116,151,169,185,205]. 

3. Early times: initial singularity 

An important question is how the higher-dimensional bulk effects modify the picture of gravitational collapse/ 
singularities. The generalized Raychaudhuri equation governs gravitational collapse and the initial singularity on 
the brane. The local energy density and pressure corrections, j^n 4 p(2p + 3p), further enhance the tendency to 
collapse [165] if 2p + 3p > (which is satisfied in thermal collapse). The nonlocal term, which is proportional to U, 
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can act either way depending on its sign; the effect of a negative U is to counteract gravitational collapse, and hence 
the singularity can be avoided in this case. 

A unique feature of brane cosmology is that the effective equation of state at high densities can become ultra stiff. 
Consequently matter can overwhelm shear for equations of state which are stiffer than dust, leading to quasi- isotropic 
early expansion of the universe. For example, in Bianchi I models [37,38,164] the approach to the initial singularity 
is matter dominated and not shear dominated, due to the predominance of the matter term p 2 /2A 2 relative to the 
shear term So/° 6 ( an d & 2 /H 2 — > as t — > 0). The fact that the density effectively grows faster than 1/a 6 for 7 > 1 
is a uniquely brane effect (i.e., it is not possible in GR). 

Indeed, a spatially homogeneous and isotropic non-gencral-relativistic brane-world (without brane tension; U = 0) 
solution, first discussed by Binetruy, Deffayet and Langlois [27,28,83], is always a source/repeller for 7 > 1 in the 
presence of non-zero shear. This solution is self-similar, and is sometimes referred to as the Brane-Robertson- Walker 
model (BRW) [51]. The equilibrium point corresponding to the BRW model is often denoted by Tb- The BRW model, 
in which 

a(t)~t^, (103) 

is valid at very high energies (p 3> A) as the initial singularity is approached (t — ► 0). In the brane-world scenario 
anisotropy dominates only for 7 < 1 (whereas in GR it dominates for 7 < 2), and therefore for all physically relevant 
values of 7 the singularity is isotropic. 

If p tot > and p tot + 3p tot > for all t < to, then from the generalized Raychaudhuri equation (51) and using the 
generalized Friedmann equation (60) and the conservation equations, it follows that for do > (where ao = a(io)) 
there exists a time tb with tb < to such that a(tb) = 0, and there exists a singularity at tb, where we can rescale time 
so that tb — and the singularity occurs at the origin. We can find the precise constraints on A and Uo in terms 
of ao in order for these conditions to be satisfied at t = to- It follows that if these conditions are satisfied at t = to 
they are satisfied for all < t < to, and a singularity necessarily results. These conditions are indeed satisfied for 
regular matter undergoing thermal collapse in which the the local energy density and pressure satisfy p(2p + 3p) > 
(and is certainly satisfied for perfect fluid matter satisfying the weak energy condition p > and a linear barotropic 
equation of state with 7 > 1). On the other hand, it is known that a large positive cosmological constant A or a 
significant negative nonlocal term hi counteracts gravitational collapse and can lead to the singularity being avoided 
in exceptional circumstances. 

In more detail, we shall see later that since the dimensionless variables are bounded close to the singularity it follows 
from Eq. (175) that < q < 2, where q is the deceleration parameter. Hence, from Eq. (176), H diverges as the 
initial singularity is approached. At an equilibrium point q = q* , where q* is a constant with < q* < 2, so that from 
Eq. (176) we have that H — > (1 + g*) -1 * -1 as t — > 0+ (r — > -00). From Eqs. (175), (177) and the conservation laws 
it then follows that p — > 00 as t — > + . It then follows directly from the conservation laws (180) that fib = K 2 p 2 /6XH 2 
dominates as t — > + (and that all of the other are negligable dynamically as the singularity is approached). The 
fact that the effective equation of state at high densities become ultra stiff, so that the matter can dominate the 
shear dynamically, is a unique feature of brane cosmology. Hence, close to the singularity the matter contribution is 
effectively given by p tot — p 2 /2X = pb, p tot = (p 2 + 2pp)/2X = (2j — l)jlb- Therefore, as the initial singularity is 
approached, the model is approximated by an RW model in some appropriately defined mathematical sense. Goodc 
ct al. [100,101] introduced the rigorous mathematical concept of an isotropic singularity into cosmology in order 
to define precisely the notion of a "Fricdmann-like" singularity. Although a number of perfect fluid cosmologies are 
known to admit an isotropic singularity, in GR a cosmological model will not admit an isotropic singularity in general. 
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V. BIANCHI IX BRANE- WORLD COSMOLOGIES 



The brane energy-momentum tensor, including both a perfect fluid and a minimally coupled scalar field, is given 



by 



where 



T — TPf -J- T s f 
± ft" — J- T -L fin , 



~ 9nv + V{<t>) 



and is the fluid 4- velocity and <j> is the minimally coupled scalar field having potential V(<j>). If 
then a scalar field with potential V(4>) is equivalent to a perfect fluid having an energy density and pressure 



(104) 

(105) 
(106) 

is timclikc, 



sf 



n s f 



(107) 
(108) 



The local matter corrections to the Einstein equations on the brane are given by Eq. (3), which is equivalent 



to 



for a perfect fluid and 



qpf 



1 2 1 . 

— p u^u u + — p (p + 2p) , 



(109) 



osf = i 

6 



1 

+ 12 



(110) 



for a minimally coupled scalar field. If we have both a perfect fluid and a scalar field and we assume that the gradient 
of the scalar field, </> ;M , is aligned with the fluid 4-velocity, u 11 (that is, </> ; ^ j \J— <f>- a & a — the local brane effects 
due to a combination of a perfect fluid and a scalar field are then 



s -1 

^ _ 12 



1 



+ v(4>) 
> a + v 



P + 2p--4>, a ^ a ~v{4>)) V- 



(in) 



The non-local effects from the free gravitational field in the bulk are given by [161,162]. We will assume 
here that D^U = = V^ v = 0. Since V^v = 0, in this case the evolution of £^ v is fully determined. In general 
U = U(t) 7^ (and can be negative) [27,28,83]. All of the bulk corrections mentioned above may be consolidated into 
an effective total energy density and pressure as follows. The modified Einstein equations take the standard Einstein 
form with a redefined energy- momentum tensor (assuming A = 0) according to Eqs. (4) and (5), where the redefined 
total energy density and pressure due to both a perfect fluid and a scalar field are given by 



„total 



total 



P + 



£ 4 



K 

12 



V(<f>) 

-1*. 

V(<f>) 



+ v(4>)) +u 



(112) 



12 ( P - o </>; <*& a + V(4>) )[p + 2p--^ a 0> a - V{<j>) 
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3 U 



(113) 



21 



As a consequence of the form of the bulk energy- momentum tensor and of Z 2 symmetry, it follows [208,212] that 
the brane energy-momentum tensor separately satisfies the conservation equations (where we assume that the scalar 
field and the matter are non- interacting) ; i.e., 

p + 3H(p + p) = 0, (114) 
• dV 

</> + 3H0+— =0, (115) 

whence the Bianchi identities on the brane imply that the projected Weyl tensor obeys the constraint 

U + 4HU = 0. (116) 



A. Bianchi IX Models: Setting up the Dynamical System 

The source term (restricted to the brane) is a non-interacting mixture of non-tilting perfect fluid ordinary matter 
with p = (7 — l)p, and a minimally coupled scalar field with an exponential potential of the form V((f>) = Voe Kk ^ 
(where (f> = <fi(t)) [103,139,206,213]. The variables are the same as those introduced in [221], with the addition of 

~ K 2 p - A ~ _ k 2 V 

3D 2 ' A ~3D 2 ' 3D 2 ' 

; |3k^ - U /-,-,„n 

*=V23D' ^ = 3D 2 

where 



4 

The total equivalent dimensionless energy density due to all sources and bulk corrections is 

2 total 2 . . 2 

ftotal = P = o + n A + $ + ^> 2 + c 2 Q u + -D 2 (O + $ + § 2 j , (118) 

where c 2 = k 4 /k 4 . 

The governing differential equations for X = [D, H, Ei, E 2 , Ai, A 2 , A3, f2, £7a, ®, &u] are as follows [113]: 



D' = -(l + q)HD 


(119) 


H' = q(H 2 - 1) 


(120) 


X' + = (q-2)H£+-S+ 


(121) 


E'_ = (q - 2)ffE_ - S*_ 


(122) 


N[ = N 1 (Hq-A±+) 


(123) 


N' 2 = N 2 {Hq + 2t+ + 2V3S_) 


(124) 


A3 = N 3 (Hq + 2E+ - 2\/3E_) 


(125) 


Q' = M(2(q+1) - 3 7 ) 


(126) 


0^ = 2HCl A (q + l) 


(127) 






= (q-2)M- ^-M> 


(128) 


$' = 2$ ^(1 + q)H + ^Y k ^j 


(129) 


£l' u = 2M u (q-l) 


(130) 



where E + , E_ are the shear variables, Ni, N 2 , A3 are the curvature variables (relative to a group- invariant orthonormal 
frame), and a logarithmic (dimensionless) time variable, r, has been defined. The quantity q is the deceleration 
parameter and S + and S- are curvature terms that are defined by the following expressions: 
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2t\ + 2t 2 _ + (37 2) n - 



+c 2 h u + -D 2 



A + 2* 2 - $ 



(ft + * 2 + l>) ((3 7 - 1)0 + 5* 2 - 



5+ 



= \ (N 2 - - i JVi (2 Aq - A> 2 - A> 3 ) 
= J V3 (iV 2 - A> 3 ) (-ATj + N 2 + A> 3 ) 



In addition, there arc two constraint equations that must also be satisfied: 

d(X) = H 2 + ^(NiN 2 N 3 ) 2 ^ 3 -1 = 



G 2 (X) = 1 - t 2 



total 



v = 



(131) 

(132) 
(133) 

(134) 
(135) 



where 



V = T2 



(N 2 + N 2 + N 2 - 2(N!N 2 + A\A> 3 + N 2 N 3 ) + 3(A>iA> 2 A> 3 ) 2 / 3 



Equation (134) follows from the definition of D, and equation (135) is the generalized Friedmann equation. The 
resulting Bianchi type IX brane- world equations are now closed and suitable for a qualitative analysis using techniques 
from dynamical systems theory. The system of equations X' = F(X) are subject to the two constraint equations 
Gi(X) = and G 2 (X) = 0. These constraint equations essentially restrict the dynamics of the dynamical system 
X' = F(X) to lower dimensional surfaces in the state space. In principle, these constraint equations may be used to 
eliminate two of the twelve variables provided the constraints are not singular. 



1. Comments 

Since the dynamical system is invariant under the transformation $ — > — $ we can restrict our state space to 
D > 0. Note that the dynamical system is also invariant under the transformation S_, Aq, N 2 , N3) — » (— — 

— N 2 , N 3 , Ni). This symmetry implies that any equilibrium point with a non-zero E± term, will 
have two equivalent copies of that point located at positions that are rotated through an angle of 27r/3 and centered 
along a different axis of the N a . 

If we assume the weak energy condition for a perfect fluid (i.e., p > 0), then we restrict the state space S to the 
set n > 0. Since we are investigating the behaviour of the Bianchi type IX brane-world models, we can restrict the 
state space to N a > without loss of generality. There are six matter invariant sets in addition to the invariant sets 
associated with the geometry of the spacetime [113]. The evolution equation (130) for Clu implies that the surface 
flu = divides the state space into three distinct regions, U + = {X e S\tlu > 0}, 14° — {X e S\&u = 0}, and 

u- = {xe s\n u < 0}. 

From (134) we have that — 1 < H < 1 and NiN 2 N 3 < 8. Furthermore, in the invariant sets U + and U°, using 
equation (225), it can be shown that 

< t 2 +1 t 2 _,n,n Al y 2 ,$,v < 1. 

However, knowing that < V < 1 and < NiN 2 N 3 < 8 is not sufficient to place any bounds on the -/V Q 's or D. 
Furthermore, in the invariant set U~ we cannot place upper bounds on any of the variables without some redefinition 
of the dimensionless variables (117). It is possible to show that the function 

W = (H 2 — 1)[- 2 "- 3 ^] fj{j»+(37-2)/3] ^ 2/ 3 (136) 

(where a, (3 are arbitrary parameters) is a first integral of the dynamical system; that is, W = [113]. In the invariant 
set U + UU°, it is possible to show that q > — 1, which implies that J7a — > as r — > — 00 for those models that expand 
for all time. Using invariants, we find that H 2 — > 1 as t — > 00, so that for ever-expanding models — > when 
7 > 4/3. We also note that the function defined by 

Z = (N^Ns) 2 , (137) 
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satisfies the evolution equation Z' = 6qZ, and is consequently monotone close to the singularity [221]. 

As mentioned earlier, generically the Bianchi type IX models have a cosmological initial singularity in which p — > oo, 
and consequently ilj, dominates as t — > + . This can be proven by more rigorous methods [196,201,202]. It can also 
be shown by qualitative methods that the spatial 3-curvature is negligable at the initial singularity. 



B. Initial Singularity 

We include a general potential for the scalar field, V, and study what happens as D — > oo at the initial singularity. 
We also include normal matter with 1 < 7 < 2. Equations (128) and (129) become 

*' = (q - 2)M - e$ (138) 
$' = 2(g+l)M + 2el>§ (139) 

where e (related to the usual inflationary slow roll parameter e) is defined by e = 3V^/2kV. 
From the Friedmann equation we have that 

ft A = — D 2 (ft + $ + * 2 ) 2 < 1 (140) 

and hence each term Dfl, DQ, D^> 2 is bounded (since the left-hand side is the sum of positive definite terms). Hence, 
as D — * 00, f2, <5 2 , $ — > 0. It is easy to show that Qa, &u - » as D — > 00. Hence (131) becomes 

q = 2{t\ + t 2 _)+AD 2 , (141) 

where 

r 2 

A= — [0 + * 2 + $][(3 7 - l)ft + 5* 2 -$] (142) 

Assuming H > 0, equations (119) and (120) imply that as r — > —00 for D — > 00, cither g — > or g is positive in 
a neighbourhood of the singularity (g can oscillate around zero; indeed, it is the possible oscillatory nature of the 
variables that causes potential problems). However, if q — ► 0, Eq. (126) implies 

0' = HVt{2 - 37) (143) 

which implies a contradiction for 7 > | (i.e., Vl -f* as r — > —00). Hence, as r — > —00, g > 0, where £) — > 00 and 

H' = -q(l-H 2 ) (144) 

and hence if — > 1 (assuming positive expansion) monotonically. [Note that this implies the existence of a monotonic 
function, and hence there are no periodic orbits close to the singularity near the set H = 1 — all orbits approach 
H = 1.] In addition, Eq. (134) gives (N^^) -> 0. 

From Eqs. (126), (127), (128) and (129) we have that as r -» -00 

I - 0, (145) 
that is, the scalar field becomes effectively massless, and 

P^0 if 7 <2. (146) 

This follows directly from the evolution equation in the case of an exponential scalar field potential, e = ^/ffc, and 
follows for any physical potential for which e is bounded as r — > — 00. Hence we obtain 

<? = 2(S 2 h + S 2 _) + (3T- 1)C 2 (147) 
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(where T = 2, C 2 = %L> 2 * 4 if 7 < 2; T = 2, C 2 = %£> 2 (0 + * 2 ) 2 if 7 = 2; T = 7 and C 2 = ^D 2 & if there is no 
scalar field), where Q\ — ► C 2 as t — > —00, and the Friedmann equation becomes 

1 - (£ 2 + S 2 .) - C 2 = -^((^1 + JV? + ) - 2(^Vi^2 + NiN 3 + N 2 N 3 )) > (148) 
From Eqs. (126) and (128) we then obtain 

(C 2 )' = 2C 2 i?[g-(3r-l)]. (149) 

We still have the possibility of S± or Ni oscillating as the singularity is approached (as in the Mixmaster models). 
A rigorous proof that oscillatory behaviour does not occur can be presented using the techniques of Rendall and 
Ringstrom [196,201] (using analytic approximations to the Brane-Einstein equations for £±, N a ; i.e., estimates for 
these quantities that hold uniformly in an open neighbourhood of the initial singularity). We can then prove that 
S ± — > as t — > —00, Q\ — ► 1, and we obtain the BRW source. 

To determine the dynamical behaviour close to the initial singularity, we need to examine what happens as D — > 00. 
Let us present a heuristic analysis (and include ordinary matter). We define a new bounded variable 

d=^j, 0<d<l (150) 

and we examine what happens as d — > 1 (assuming H > 0) . From above we have that d — ► 1 and H — ► 1 monotonically, 
and hence we need to consider the equilibrium points in the set d = 1. 

The analysis depends on whether the quantity A defined by (142) that occurs in the expression for q in equation 
(141) is zero or not in an open neighbourhood of the singularity. Without loss of generality, assuming that A > 0, we 
define a new time variable by 

> = (!-<*)% ; 
AH 

and the remaining evolution Eqs. (on d = 1) become 



£+ = £+, Ni=N u N 2 =N 2 , N 3 = N 3 , 

n = 20, = 20 A , * = $ = 2$, n u = 2n u . (151) 

Therefore, the only equilibrium point is 

S+ = E_ = TVi = N 2 = N 3 = n = O a = * = $ = = 
and this is a local source. This equilibrium point corresponds to the BRW solution [50] with 

r 2 

— L> 2 * 2 = 1, D 2 $ = (L> 2 = 0). 

In the next two sections we will consider the two cases of a scalar field and a perfect fluid separately. 

Due to the quadratic nature of the brane corrections to the energy momentum tensor, a rich variety of intermediate 
behaviour is possible in these two fluid models. We note that a dark radiation density U together with a scalar field 
is similar to previous scaling models [26,114]. Here the equivalent equation of state is pu = \U. It is known that the 
bifurcation value for these scaling models is k 2 — 37, which for 7 = 4/3 corresponds to a value k 2 = 4, a bifurcation 
value found in the analysis here. The intermediate behaviour of these multifluid models is extremely complex [113]. 
A more complete analysis of the Bianchi type II models (in which some of the intermediate behaviour is outlined) 
was studied in [109,111]. 

VI. QUALITATIVE ANALYSIS 
A. Scalar Field Case 

In the invariant set U + UU°, we can show that q > — 1. This implies that the set ft\ = is the invariant set 
containing all of the past asymptotic behaviour for ever-expanding models in U + UU°. It can be argued that a scalar 
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field becomes essentially massless as it evolves backwards in time, hence it will dominate the dynamics at early times. 
In an effort to understand the dynamical behaviour at early times we shall first assume that there is no perfect fluid 
(O = 0) and the 4D cosmological constant is zero. 

If f2 = Oa = (i.e., for a scalar field only), then the dynamical system simplifies (effectively the evolution equations 
for Cl and Q,\ are omitted). The curvature terms 5+ and S- arc defined as before, and the deceleration parameter q 
becomes 



q = 2E 2 + 2Y?_ + 2* 2 - $ + c 2 



^D 2 (V 2 + $>)(5^ 2 - + Q u 



(152) 



The two constraint equations are given explicitly by 



1 = H 2 + -(N!N 2 N 3 ) 2/3 



H 2 



£ 2 + + £ 2 _ + h total + ((iV 2 + N 2 + N 2 ) - 2(^2 + N 2 N 3 + 7ViiY 3 )) 



(153) 
(154) 



where 



Qtotal _ K P 



2 Jtotal 



3D 2 



= $ + f 2 + c 2 



D 2 (<s> + v 2 y + n u 



(155) 



We first consider the equilibrium points at finite D. We define X = {D, H, E+, E_, iVi, N2, N$, *, $, Qu} and we 
restrict the state space accordingly to be S — {X e 5|0 = 0, 17a = 0}. (Note, e is a discrete parameter where e = 1 
corresponds to expanding models, while e = — 1 corresponds to contracting models.) There are a number of saddle 
points [113]. The zero curvature, power law inflationary Robertson- Walker equilibrium point P f , defined by 

X = [0, e, 0, 0, 0, 0, 0, — ^p^, 1 — 0] is a local sink in the eight dimensional phase space S if k 2 < 2 (if e = — 1, this 
point is a local source). When 2 < k 2 this point becomes a saddle. We note that there are models that recollapse. 
There are no sources for finite values of D. 

We considered the equilibrium points at infinity in the previous section. We now show that (for no perfect fluid 
and a scalar field with an exponential potential) the BRW solution is always an equilibrium point of the system and a 
local stability analysis establishes that it is a local source. In order to analyze the dynamical system for large values 
of I?, we define the following new variables 

$ = r 2 sin 2 0, * = rcos0, Cl x - i C 2 D 2 ($ + * 2 ) 2 = Jc 2 D 2 r 4 , 

so that the variable D is essentially replaced by the bounded variable £7a in the set UUU° (hence in this set the only 
variables that remain unbounded are the Net's). The dynamical system becomes 



& = q(H 2 - 1) 


(156) 


S' + = (q - 2)Ht + - S+ 


(157) 


S'_ = {q- 2)#S_ - S- 


(158) 


N[ = N 1 (Hq-At+) 


(159) 


N' 2 = N 2 (Hq + 2S+ + 2V3S_) 


(160) 


= N 3 (Hq + 2S+ - 2\/3S_) 


(161) 


% = 2m u {q-l) 


(162) 


Q! x = 2HCl x (q+l - 6 cos 2 6>) 


(163) 


r' = rH(q- 2 + 3sin 2 6>) 


(164) 


6' = sind(3H cos 9 +^-kr) 


(165) 



where 



q = 2Y? + + 2Y?_ + r 2 (3 cos 2 6 - 1) + Cl x (6 cos 2 6-1) + c 2 £l u 



(166) 
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The two constraint equations become 

l = ff 2 + i(A>iA 2 A 3 ) 2/3 (167) 

h 2 = t\ + 1 2 _ + n x + c 2 n u + r 2 

+ T2 + ^ 2 + ^ ~ 2 ^^ 2 + ^ 3 + ^^ 3 )) ( 168 ) 

From the constraint equations, we have < {H 2 , S+, E 2 ,, fiw, ^a, r} < 1, < 6 < ir, and N1N2N3 < 8 (assuming 
n u > 0). 

The BRW solution is represented by an equilibrium point in the set ft\ ^ (D — > 00). If we let X = 
{H, E + , E_, Ai, A2, A 3 , £7^, r, then the relevant equilibrium points are X = [e, 0, 0, 0, 0, 0, 0, 1, 0, |± §]. Using 
the constraint equation to eliminate f2y , the eigenvalues of the linearization at the points 9 — 0, it are 

e(10,10,5,5,5,3,3,3,3). 

A value of 9 that satisfies 9' = in a neighbourhood of an equilibrium point corresponds to a tangent plane to an 
invariant surface passing through that equilibrium point. In the analysis above, the directions 9 — and 9 = tt 
correspond to the $ = invariant surface (i.e., the massless scalar field models). We observe that this equilibrium 
point is a source that strongly repels away from ^ = (for e > 0). That is, when traversed in a time reverse direction, 
typical orbits would asymptotically approach a massless scalar field BRW solution. 

In summary, the isotropic BRW solution is a global source and the past asymptotic behaviour of the Bianchi IX 
brane world containing a scalar field is qualitatively different from what is found in GR. 

We note that assuming Qu > 0, the future asymptotic behaviour of the Bianchi IX brane world containing a scalar 
field with an exponential potential is not significantly different from what is found in GR [113]. For < k < \/2, 
the future asymptotic state of ever-expanding models is characterized by the power-law inflationary solution, and if 
k > \/2 there no longer exists any equilibrium point representing an expanding model that is stable to the future. We 
therefore conclude that if k > \/2, then the Bianchi IX models must recollapse. In [112] it was shown that if k > \/2, 
then a collapsing massless scalar field solution is a stable equilibrium point. In the brane- world scenario, we have that 
this final end-point is the BRW solution. However, if VLu < 0, then a variety of new behaviours are possible, including 
possible oscillating cosmologies [38]. 

B. Perfect Fluid Case: No Chaos in Brane- World Cosmology 

We shall now discuss the asymptotic dynamical evolution of spatially homogeneous brane- world cosmological models 
with no scalar field close to the initial singularity in more detail. Due to the existence of monotone functions, it is 
known that there are no periodic or recurrent orbits in orthogonal spatially homogeneous Bianchi type IX models in 
GR. In particular, there are no sources or sinks and gcnerically Bianchi type IX models have an oscillatory behaviour 
with chaotic-like characteristics, with the matter density becoming dynamically negligible as one follows the evolution 
into the past towards the initial singularity. Using qualitative techniques, Ma and Wainwright [221] have shown that 
the orbits of the associated cosmological dynamical system are negatively asymptotic to a lower two-dimensional 
attractor. This is the union of three ellipsoids in R 5 consisting of the Kasncr ring joined by Taub separatrices; the 
orbits spend most of the time near the self-similar Kasner vacuum equilibrium points. More rigorous global results are 
possible. Ringstrom has proven that a curvature invariant is unbounded in the incomplete directions of inextendiblc 
null geodesies for generic vacuum Bianchi models, and has rigorously shown that the Mixmaster attractor is the past 
attractor of Bianchi type IX models with an orthogonal perfect fluid [201,202]. 

All spatially homogeneous Bianchi models in GR expand indefinitely except for the type IX models. Bianchi type 
IX models obey the "closed universe recollapse" conjecture [155,156], whereby initially expanding models enter a 
contracting phase and recollapse to a future "Big Crunch" . All orbits in the Bianchi IX invariant sets are positively 
departing; in order to analyse the future asymptotic states of such models it is necessary to compactify phase-space. 
The description of these models in terms of conventional Hubble- or expansion-normalized variables is only valid up 
to the point of maximum expansion (where H = 0). An appropriate set of alternative normalised variables leading to 
the compactification of Bianchi IX state space, which were suggested in [221] and utilized in [53,62,112,131,218], was 
discussed in the previous section. 

Let us discuss the Bianchi type IX brane-world models with a perfect fluid source. Unlike earlier, and following 
[37,38,221], we define Hubblc-normalized shear variables E + , E_, curvature variables Ni 7 N 2 , A3 and matter variables 
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(relative to a group-invariant orthonormal frame), and a logarithmic (dimensionless) time variable, t, defined by 
dr = Hdt. These variables do not lead to a global compact phase space, but they are bounded close to the singularity 
[221]. The governing evolution equations for these quantities are then 

£'+ = (q- 2)S+ - S+ (169) 

£'_ = (g - 2)£_ - S_ (170) 

N[ = (q- 4E+)iVi (171) 

N' 2 = (g + 2£+ + 2\/3£_)jV2 (172) 

N' 3 = (q + 2S + -2a/3£_)7V 3 , (173) 

where a prime denotes differentiation with respect to r, and K, S+ and S- arc curvature terms that are defined as 
follows: 

K = h. ^ + N * + N ^ ~ 2 ^ NiN2 + N2Ns + NiN ^) ( l74a ) 

S+ = i (N 2 - N 3 ) 2 - i Ni (2 Nx - N 2 - N 3 ) (174b) 

S- = i V3 (N 2 - N 3 ) (-JVi + iV 2 + jV 3 ) (174c) 
The quantity g is the deceleration parameter given by 

g = 2S^ + 2S 2 _ + i^^ + ^P,. (175) 

The decoupled Raychaudhuri equation becomes 

H' = -(l + q)H. (176) 

In addition, the generalized Friedmann equation reduces to the constraint 

1 = + X? + Sli + K. (177) 

Due to the symmetries in the dynamical system, we can restrict ourselves to the set N\ > 0, N 2 > 0, and N 3 > 0, 
without loss of generality. We again note that these new normalized variables, for which the evolution equation for 
H has decoupled from the remaining equations, are bounded close to the initial singularity. 
In the above 

^.2 „tot 2 ntot 

where the fti (i = 1 — 4) are given by 

2 

ni=fi = ^; P=(7-l)fi (179a) 



3H 2 

2 2 

" 2 = 0fe = 6Tff2 ; p " = ( 2 t - ( 179b ) 



07V 1 

° 3 = " W = T^Atf^ Pw = 3°" (179C) 
" 4 = " A = 3^ 5 Pa = ~" A ' (179d) 

n{ = [2(1 + g) - Srjfii, (180) 

ri = 7,r 2 = r fe = 2 7 ,r 3 = r w = |r 4 = r A = o. (isi) 

An extensive analysis of scalar field models has shown that close to the initial singularity the scalar field must be 
massless [53,62,84,197]: it is therefore plausible that scalar field models can be approximated by a stiff perfect fluid 
close to the initial singularity (particularly regarding questions of stability). 



which satisfy the equations 



where 
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1. The isotropic equilibrium point 



We recall that generically the Bianchi type IX models have a cosmological initial singularity in which p — > oo, and 
consequently fif, dominates, as t — > + . It can also be shown by qualitative methods that the spatial 3-curvature is 
negligable at the initial singularity. Hence, close to the singularity p tot = pb, so that = ^6 and ^2 Pi = Pb, 

where Pb = (2-f — l)f2&. Consequently, Eq. (177) can be written as 

n b = l- Y? + -Y?_- K, (182) 

which can be used to eliminate from the governing equations. In particular, Eq. (175) becomes 

q = 3(1 - 7 )(S 2 + + E 2 _) + (37 - 1)(1 - K), (183) 

and the governing equations are given by the dynamical system above, where q is now given by Eq. (183). 

There is an equilibrium point of the dynamical system, denoted by T bl given by S + = £_ = 0, and TVi = iV 2 = 
N 3 = 0, which corresponds to the BRW model. For the equilibrium point Tb, the 5 eigenvalues are: 

3( 7 - 1), 3( 7 - 1), (3 7 - 1), (37 - 1), (3 7 - 1) (184) 

Hence, for all physically relevant values of 7 (7 > 1), T b is a source (or past-attractor) in the brane-world scenario 
and the singularity is isotropic. This contrasts with the situation in GR in which anisotropy dominates for 7 < 2. 
This is also consistent with previous analyses of Bianchi type I and V models where T b is always a source for 7 > 1 
[37,38]. 

Therefore, in general the Bianchi IX brane-world models do not have space-like and oscillatory singularities in the 
past, and consequently brane-world cosmological models do not exhibit Mixmaster and chaotic- like behaviour close to 
the initial singularity. We expect this to be a generic feature of more general cosmological models in the brane-world 
scenario. Following arguments similar to those above, it follows that there exists a singularity of a similar nature in 
all orthogonal Bianchi brane-world models. In general, Sib will again dominate as the initial singularity is approached 
into the past and the qualitative results will follow. In particular, T b will be a local source in all spatially homogeneous 
models and in general the initial singularity will be isotropic. 



29 



VII. INHOMOGENEOUS BRANE- WORLD MODELS 



We have found that an isotropic singularity is a past-attractor in all orthogonal Bianchi models. It is plausible that 
the cosmological singularity is typically isotropic in spatially inhomogeneous models [50] . We shall study this further 
by considering the dynamics of a class of G2 spatially inhomogeneous cosmological models with one spatial degree of 
freedom in the brane- world scenario [57]. We follow the formalism of [81], which utilizes area expansion normalized 
scale-invariant dependent variables, and we use the timelikc area gauge to discuss the asymptotic evolution of the 
class of orthogonally transitive G2 cosmologies near the cosmological initial singularity. The initial singularity can 
be shown to be characterized by Ei 1 — > as t — > —00, where Ei 1 is a normalized frame variable [81]. In particular, 
we shall discuss the numerical behaviour of the local dynamical behaviour of this class of spatially inhomogeneous 
models close to the singularity. 

All of the bulk corrections may be consolidated into an effective total energy density, pressure, anisotropic stress 
and energy flux as follows. The modified Einstein equations take the standard Einstein form with a redefined energy- 
momentum tensor (possibly tilting now): 



(185) 



where 



= k z p + 



V = 



K 2 p + 



6^ 
A 
6k 2 



A 



(2p 2 - 3 + —U 



1 1 



— (2p 2 + App + V ^ - Aq^) + 



6 k 2 



A 

6k 2 r 
A 



^ (2pq M - 3n^q") + -^Q^ 



(186) 
(187) 
(188) 
(189) 



We shall assume that V^ v = 0, so that the evolution of £^ u is fully determined [208,212]. In the inhomogeneous 
cosmological models of interest here a non-zero D M p acts as a source for Q M , and hence — is not consistent with 
an inhomogeneous energy density and we need to include a dynamical analysis of the evolution of Q M . We shall make 
no further assumptions on the models and include all terms in the numerical analysis. 



2. Initial singularity 

From the numerical analysis we find that the area expansion rate increases without bound (f3 — ► 00) and the 
normalized frame variable vanishes (Ei 1 — ► 0) as logarithmic time r — > —00. Since j3 — > 00 (and hence the Hubble 
rate diverges), there always exists an initial singularity as r — > — 00. Thus the singularity is characterized by Ei 1 = 0, 
which allows both dynamical and numerical results to be obtained. As noted earlier, since p — > 00 as t — > —00, 
it then follows directly from the conservation laws that pb ~ p 2 dominates as r — > —00 and that all of the other 
contributions to the brane energy density are negligible dynamically as the singularity is approached. Hence close to 
the singularity the matter contribution is effectively given by p tot = (2-y — l)p tot . 

Models with an isotropic initial singularity [100,101] satisfy lim r ^_ oc fi^ = 1, lim T _>_ 00 v = 0, lim T ^_ oc E 2 = 0. 
Their evolution near the cosmological initial singularity is approximated by the BRW model corresponding to the 
'equilibrium point' T\>, characterized by 

n b = 1; = E x x = £+ = £_ = E x = A_ = N x = v = Q u = Q = fl u . (190) 

It was argued above that for all physically relevant values of 7 (7 > 1) J 7 ;, is a source (or past-attractor), and hence the 
singularity is isotropic, in all non-tilting spatially homogeneous brane- world models [50]. It was also shown that Ti> is a 
local source in the family of spatially inhomogeneous 'non-tilting' G 2 cosmological models for 7 > 1 in [50] . Linearizing 
the evolution equations about !Fi,, using the same spatial reparametrisation as in [81] (so that Ei 1 — exp((37— l)r)), 
the following general solution of the linearized equations was obtained [50] : 

E_=ai(x)exp(3(7-1)T), E x = a 2 {x) exp (3( 7 - l)r), (191) 
A_ = a 3 (x) cxp ((3 7 - l)r), N x = a 4 (x) exp ((3 7 - l)r), (192) 
n b = 1 + a 5 (x) exp (3(7 - 1)t), v = a 6 (x) cxp ((37 - l)r) (193) 
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where the cti(x) are arbitrary functions of the space coordinate. In addition, it follows from the conservation laws 
that Q — > 0, flu — > and — > as the initial singularity is approached. The above linearized solution represents 
a general solution in the neighbourhood of the initial singularity. Hence Tb is a local source or past-attractor in this 
family of spatially inhomogeneous cosmological models for 7 > 1. In particular, we see that the shear and curvature 
asymptote to zero as r — > —00, and hence the singularity is isotropic. We also note that, unlike the analysis of the 
perfect fluid GR models in [81], the Kasner equilibrium set K are found to be saddles in the class of G 2 brane-world 
cosmological models. 



A. G2 Brane Cosmology 

We can study the nature of the initial singularity in the spatially inhomogeneous brane cosmological models. An 
analysis of the behaviour of spatially inhomogeneous solutions to Einstein's equations near an initial singularity 
has been made in classical GR [23,24,87,118,130,222]. We consider the class of G2 cosmologies with two commuting 
Killing vector fields, which conscqucnly admit one degree of spatial freedom [81] . The evolution system of Einstein field 
equations are partial differential equations (PDE) in two independent variables. The orthonormal frame formalism is 
utilized [80,166] with the result that the governing equations are a first-order autonomous equation system, and the 
dependent variables are scale- invariant. In particular, we define scale-invariant dependent variables by normalisation 
with the area expansion rate of the G2-orbits in order to obtain the evolution equations as a well-posed system of 
PDE, ensuring the local existence, uniqueness and stability of solutions to the Cauchy initial value problem for G2 
cosmologies. Following [81] we assume that the Abelian G2 isometry group acts orthogonally transitively on spacclikc 
2-surfaces, and introduce a group- invariant orthonormal frame { e a }, with e 2 and e3 tangent to the G2~orbits. The 
frame vector field eo, which defines a future-directed timelike reference congruence, is orthogonal to the G2-orbits 
and it is hypersurface orthogonal and hence is orthogonal to a locally defined family of spacelike 3-surfaces. We then 
introduce a set of symmetry-adapted local coordinates { r, x, y, z} 

e = A -1 d T , ei = ei 1 d x , = ca b x b , A, B = 2, 3 , 

where the coefficients are functions of the independent variables t and x only. The only non-zero frame 
variables are thus given by A , ei 1 , eA B , which yield the following non-zero connection variables: 
a, [3, 01, n+, er_, n x , cr x , n,-, in, fii. The variables a, f3, o~ and a x are related to the Hubble volume expan- 
sion rate H and the shear rate a a p of the timelike reference congruence e ; in particular, 6 := 3H = a + 2/3. The 
variables ai, n+, n x and n_ describe the non-zero components of the purely spatial commutation functions a a and 
n a p [221]. Finally, the variable iii is the acceleration of the timelike reference congruence e , while VLi represents the 
rotational freedom of the spatial frame { e a } in the (e 2 , e 3 )-plane. Setting Q 1 to zero corresponds to the choice of a 
Fermi-propagated orthonormal frame { e a }. Within the present framework the dependent variables 

{ A, in, fii } (194) 

enter the evolution system as freely prescribable gauge source functions. 

Since the Gi isometry group acts orthogonally transitively, the 4-velocity vector field u of the perfect fluid is 
orthogonal to the G2-orbits, and hence has the form 

u = r(e + wei), (195) 

where T= (1-v 2 )- 1 / 2 . We assume that the ordinary matter is a perfect fluid with equation of state 

Pfl = (7 - l>/; • (196) 

In a tilted frame, we have 

= pu^u u + ph^ v + + 2q {ll u v) , (197) 

where 

G+ 1 3( 7 -l)(l-f 2 )+7^ 2 IP IP nQfiA 

p = 1 _ v2 Pfu p = 3 p, q a = ~Q~^ Va > n <*P = G^ V{aV(3) ^ ' 

and G + = 1 + (7 — l)v 2 . The basic variables are p and v a . The quadratic correction matter tensor S a p is explicitly 
given in [57]; comparing with (198), we see that it is similar to a perfect fluid with equation of state pfi = (27 — l)p/z . 
In addition, 
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£^ = — — {p u u a u u +p u h^ + ttV + 2«"(„u„)), (199) 

and is similar to a perfect fluid with equation of state p u fi = \p u fi > but with non-zero {q u fi) u and (n u ji)^- We 
shall use p u and q u M as the basic variables, and we can set (7r u /z) M1/ = self-consistently [57]. We define qt, qb, q u , 
TTt , and tt u as in [57] , so that the orthonormal frame version of the field equations, matter equations and non-local 
equations, when specialised to the orthogonally transitive Abelian G 2 case with the dependent variables presented 
above, can be written down. 



1. Scale-invariant reduced system of equations 

We introduce /3-normalised frame, connection and curvature variables as follows [57,81]: 

{M-\E l 1 ):=(N-\e l 1 )/P (200) 
{U, A, (1-3E+), £_, N x , E x , JV-, N+, R) := {u lt a u a, cr_ , n x , <r x , n_, n+, fii )//? (201) 

fin- 2 fife 2 fi/i fiO 

(fi, fi„, Q u ) := ( K 2 p, ^Pn, —qu)/(ZP 2 ) = (« 2 p, S )/(3/?2) (2 ° 2) 

s := (203) 

where (5 is the area expansion rate of the G 2 -orbits. The new dimensionless dependent variables are invariant under 
arbitrary scale transformations, and are linked to the _ff-normalised variables through the relation H = (1 — £+) /? 
[221] . Note that v is dimensionless. It is necessary to introduce the time and space rates of change of the normalisation 
factor (3; we algebraically define q and r in terms of the remaining scale-invariant dependent variables: 

Af- 1 d T B = -(q+l)B, 0= (E 1 1 d x + r)B . (204) 

The various physical quantities are defined by 

fitotai = fi + fi 6 + fi u = fi + } (1 + (2 7 - \)v 2 ) fi 2 b 2 + n u (205) 
iw = p + + p. = 3( -'-"{ G ; 2) + -'" 2 n + 3( ^-^-;;> + 2 ^ 2 (i - ,*>tf fl" + in. » 
Qtotai = Q + Qb + Qu = ^fi + 7 ^fc^ » 2 # 2 + Qu (207) 

(jr_|_ D (jt_|_ 

n„„ . n + n 4 + n, = H2^n + ^^n> b> + (208) 

The scale-invariant system contains evolution equations for the dependent variables {-Ei 1 , £+, A, N + , S_, A<x, 
S x , JV_, fif,, u, Q}, but not for the gauge source functions {AT, U, R} , which are arbitrarily prescribable real- 
valued functions of the independent variables t and x. We can fix the spatial gauge by requiring 

N + = V3N-, i?=-V3£ x , (209) 

which is preserved under evolution and under a boost. We fix the temporal gauge by adapting the evolution of the 
gauge source function U. The separable area gauge is determined by imposing the condition 

= (r - U) , (210) 

which determines U algebraically through the defining equation for r. It follows immediately from the gauge fixing 
condition that N = /(t), and we can use the r-reparametrisation to set /(r) = Mq, a constant, which we choose to 
be unity; i.e., TV = A/o := 1. Therefore, r is effectively a logarithmic proper time, and the initial singularity occurs 
for r — > — oo. From the area density A of the G2-orbits, defined by A^ 1 = e 2 2 e 3 3 — e 2 3 e^ 2 , we obtain the magnitude 
of the spacetime gradient (V a A) (V a A) = - 4/3 2 (1 - A 2 ) A 2 , so V ' a A is timelike for A 2 < 1. For the class of G 2 
cosmologies in which the spacetime gradient V a A is timelike we can choose the gauge condition A = 0, which would 
be achieved by choosing e to be parallel to V a A. It follows that d x A = 0, and we obtain A = l 2 ) c 2J ^° T , which is a 
function of r only. This is refcred to as the timelike area gauge [81]. 
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B. Governing Equations in Timelike Area Gauge 



Let us explicitly give the evolution system in the timelike area gauge. In the timelike area gauge we can use equation 
(209) to eliminate the evolution equation for N + , and the evolution equation for A becomes trivial. The relevant 
equations are: 
Evolution system: 

d T B = -(q+l)B (211) 

drE^ = ( (? + 3S+)£; 1 1 (212) 

9 T E+ = {q + 3S+ - 2) £+ - 2 (N 2 + TV 2 ) - § E^ d x r - § n total (213) 

<9 T E_ + Si 1 d x N x = (q + 3S+ - 2) £_ + 2V3£ 2 - 2V3iV 2 (214) 

d T N x +E 1 1 d x Z_ = (q + 3Z + )N x (215) 

a T S x - Ei 1 d x N_ = (q + 3S+ - 2 - 2V3S_) S x - 2^3 iV x 7V_ (216) 

d T N_ - E^ d x E x = (q + 3S+ + 2\/3S_) 7V_ + 2\/3E x iV x (217) 

^(d T + -^-vE 1 1 d x )Q + f 2 E 1 1 d x v = 2^-f 1 [^±(q + l)- i(l-3S+)(l+« 2 )- 1] (218) 

Ji G + G + 7 Z 

; r^(^- ?7 4^^^i 1 ^)«+ f 2 E 1 1 d x n = --^n(i-v 2 )[¥—^G +r 

h G + G_ /lCr- 7 



d T Q u + E^dxQy 



+ (2-7)(l-3E + )w-2( 7 -l)t;] (219) 
0, + (220) 



1 + w 2 



2(q - 1) + 12— -£+ 

3 + ir 



where 



and 



2d - ii 2 ") 

d T Qu + E^d^Pu + n„) = 2(g + 3S+ - 1)Q U - 1 2 ' rSl u + X (221) 

O i 1) 



1 _l Q ?; 2 i a v 

d x (P u + n M ) = -—^ d x n u + o M d x v , (222) 

3 + v z (3 + w^) 2 



v _ 7(1 - v 2 )B 2 n 2 
~x = ^ 

6G+ 2 G_ 



7u(l - v 2 )(l - 3S+) + 2 7 w + 2(1 - v 2 )r 



6G+ 2 G_ 6G+ 3 G_ V ' 

Constraint equations: 

= 2S+ - 1 + y?_ + s 2 + iv 2 + iv 2 + n + ^~jf} (1 + (27 - 1) « 2 ) ^ 2 £ 2 + n u (224) 

12 G+ 



jv 7 v (1 — v 2 ) 

-jv_z, x ; + 2l 

Defining equations for q and r: 



= r + 3 (7V X - 7V_ E x ) + f ( ±-Sl + ' ^ > B 2 + (.)„ ) (225) 



g :=i + |(S 2 _+iV 2 +S 2 +7V 2 ) 



Q + 3 £ + — 2 + |( — fltotal + P total + n to tal) 

,,,,, 2 fh-W^v^W n u (2- 7 )n \ 

= 2 + 3(1 - W H 12G? sT^^G^J (22?) 

r := , (228) 
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where 

G ± := 1 ± (7 - 1) v 2 , h := ^zil (1 - v 2 ) 2 , f 2 := (1 - v 2 ) 2 , f 3 := (3 7 - 4) - (7 - 1) (4 - 7 ) « 2 • (229) 



1. Numerical results 

We have written the governing equations as a system of evolution equations subject to the constraint equations 
(224) and (225). We can use (224) to obtain fl u and (225) to solve for Q u , and thus treat the governing system 
as a system of evolution equations without constraints. From the numerical analysis [57] we find that the area 
expansion rate increases without bound (fi oo) and the normalized frame variable vanishes {Ex 1 — > 0) as r — > — oo. 
Since (3 — > oo (and hence the Hubble rate diverges), there always exists an initial singularity. In addition, we find 
that {fl, E_, A^ X ,S X , N_,r} — > as t — > — oo for all 7 > 1. In the case 7 > 4/3, the numerics indicate that 
{v, fl u , Q u , E+} — > (and fib — » 1) for a/i initial conditions. In the case of radiation (7 = 4/3), the models still 
isotropize as r — > —00, albeit slowly. This degenerate case was investigated in more detail [57]. For 7 < 4/3, 
{v, fl u , Qu, 5]+} tend to constant but non-zero values as r — > —00. It is interesting to note that v 2 -/* 1; i.e., the tilt 
does not tend to an extreme value. 

The numerical results support the fact that all of the models have an initial singularity and that for the range of 
values of the equation of state parameter 7 > 4/3 the singularity is isotropic. Indeed, the singularity is isotropic for all 
initial conditions (and not just for models close to Fb), indicating that for 7 > 4/3 this is the global behaviour. The 
numerical results also support the exponential decay (to the past) of the anisotropies. Indeed, we found no evidence 
that models do not isotropize to the past for 7 > 4/3. 

In the case 7 < 4/3 we find numerically that — > 0, so that there is always an initial singularity as r — > —00. 
In addition, we find that Ei 1 , r, JV_, iV x , E_ E x , fl all vanish as the initial singularity is approached. However, the 
initial singularity is not, in general, isotropic. 

We can derive the asymptotic dynamics of the models. We have that lim T _ > _ (X) (.Ei 1 , B^ 1 , r, N-, N x , E_, S x , E+, 
fl, v, fl u , Q u , fib ~ 1) = from the numerics, and that d x {E\ 1 , B~ x , r, N^, N x , E_, E x , E + , fl, v, fl u , Q u , fi& — 1) are 
bounded as r — > —00, and V = C(/(r)), which implies that d x V = 0(f(r)), when 7 > 4/3. In particular, since 
Ei 1 — > as t — ► —00, we can follow the analysis in [154] and use equations (211)-(228) to obtain the following 
asymptotic decay rates [57]: 

E x x = 0( e K 3 r- 1 >- e ] T ) ) N x = 0(e^-^-^ T ), B- 1 =0(c^-^ T ), (230) 
=0(c^- 1) - £ ^ T ), E x =0(c^-^-^ T ), E_ = ^(el 3 ^- 1 )- 6 ^), (231) 
fl = 0{c^-^ T ), r = 0(e^- 1 ^ T ), (232) 

V = 0(C^-^ T ), E+ = 0(e^-V-^ T + 0^(37-4)-^^ Qu = 0{c K^-i)-e]r^ ^ = 0(e [2(3 7 -4)- e ]r j ^ (333) 

for any e > 0. The asymptotic dynamics for 7 = 4/3 are also given in [57]. 

C. Discussion 

1. The invariant set Ei 1 — 

From numerical experiment we find that as r — > —00, Ei 1 — > very rapidly, E_, ExjiVx,^-, r — > and 
B — > 00. Let us discuss the dynamics in the invariant set Si 1 = at early times. In order to find the early time 
behaviour of fib in the invariant set Ei 1 — 0, we obtain the evolution equation for fl b : 

d T flb = d v flb d T v + dnflb d T fl + dsflb d T B, (234) 

where d T B, d T fl and d T v are given separately by (211), (218) and (219) with Ei 1 = 0. 
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Equations (211)-(221) imply that X_ = S x = iV x = N- = fi = is an invariant subset of the invariant set Ei 1 = 0. 
In this invariant subset 

~(Q b + Qu) = r = -^^-=0 (235) 

The dynamics in the invariant subset f?:i? 1 1 = E_=X x = iV x = iV_ = fi = (in which r = Q b + Q u = 0) is 
given by 

_ n . r 2 + en b + en u - 37(1 + n b + o„) \ 24 

+ 2^ b {(, + !)- 3(1 - E + ) (1 + 3(7-l)(l-^) + 7^ } + s+ |^| (236) 
9 tW = _ v ( 1 -^ 2 ) ((2 - 7 )(1 - 3E+) - 2( 7 - 1)) (237) 

d T n u = 2( q + 1) (fi 6 + n u ) - 3(1 - e+) |fi 6 + fi u + i ( 3(27 ( ^ ) ( ( 2 1 ~_' ;2 1 ) ) ^) 7 ' ;2 ) n * + ^ 

[4 7^ 2 St; 2 \ f 2 + 6fi fc + 6fi M - 3 7 (1 + »«, + » M ) ) 2 4q 
+3S+ \ 3 (1 + (2 7 - I), 2 ) n " + 3(3T^) °« / " I (l + (2 7 -l), 2 )G_G+ } 7 " ° b 

-2 *{(* + !)- 3(1 - E + ) (1 + 3(7 - 1)( ^; 2)+7 ^ ) + (238) 

£+ = i (1 - fi b - fi„) (239) 

+ 1 f 3(27-l)(l-^) + 27^ 2 , Q 4 7 ^ 2 8^ 2 1 

9 - 2 + 2 \ (l + (2 7 -l)« 2 ) n " + ^ + (l + (2 7 -l), 2 )° b + (3T^)""/ (240) 

and we note that Q u = —Q b = —2jvfl b /(l + (27 — 1) v 2 ). 

Since fit <~ fi 2 B 2 , in this invariant set we have a closed system of ODE (236)- (238) in terms of the three variables 
fib, v and fi M . The equilibrium points are summarized in Table II in [57]. The sources are T b (fit = 1, v = 0, fi u = 0) 
for 7 > 4/3 and V± (fib = -Fi( 7 ), v = ±F 2 ( 7 ), fi„ = ^3(7)) for 7 < 4/3 (where the constants Fi, F 2 , -F3, which depend 
on 7 , are given explicitly in [57]). 

Therefore, in the invariant subset E, T b is the global source for 7 > 4/3, and V± (± depending on sign of v) are 
anisotropic sources for 7 < 4/3. There is a bifurcation at 7 = 4/3, where V± and T b coincide; from the dynamics in 
the center manifold it was found [57] that the nonhyperbolic equilibrium point is a source in the 3D invariant subset. 
This is consistent with the decay rates given above and with the numerical analysis. Hence models isotropize 'slowly' 
in the radiation case. 

In [57] it was shown that when 7 > 4/3, T b is also a source for the full state space and that when 7 < 4/3, V± arc 
sources for the full state space. We observe that at T b 



where 



q + 1 - 3 7 > 0, q + 3S+ - 2 = 3(7 - 1) > 0, Ti = I > 0, T 2 = 3 7 - 1 > 0, 



and at V± 



27 2(7-1) (37-2)7 2 

q+l = —^>0, < ? + 3 £ + -2 = ^ 7 -^ >0, T 1 = - l 2 7 n >0, T 2 = - >0, 

2 — 7 (2-7) (2 7 2 -77 + 4) 2 — 7 

where T\ and T 2 are given explicitly in [57]. So near the equilibrium points T b and V±, q+ 1 > and g + 3S+ — 2 > 
and thus as r — > —00, we have B — > +00 and i^i 1 — > 0. Near J^, and Pi, we can neglect the terms with l E\ d x ' in the 
equations (211)-(220) and treat the PDE as a system of ODE. Then from the linearization of equations (214)-(218) 
near T b and V±, and using the fact that q + 3X+ — 2 > 0, T\ > near T b and T±, we find that as t — > —00, S x , 
A^x , iV_ and fi (and hence Q) will decrease monotonically towards zero. We also observe that near T b and V± 

d T (Qb + Qu) = 2(q + 3S+ - 1) (Qb + Qu) (2(q + 3S+ - 1) - I (-fib - fi„ + P b + P u + U b + n u )j (Qb + Qu) 

+ f3(iV_£ x -7V X E_)-^Q) (-fib-fi u + Pb + P u + nb + n„). (241) 
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Linearizing this evolution equation for Qb + Q u near Tb and V± and using T 2 > 0, we conclude that Qb + Q u , and 
thus r, decrease towards zero as r — > —00. Therefore, as r — > —00, all orbits in the full state space evolve back 
towards the invariant subset E. Because Tb (when 7 > 4/3) and V± (when 7 < 4/3) are sources in this invariant 
subset, orbits will shadow orbits in the invariant subset and evolve back towards Tb (when 7 > 4/3) or T± (when 
7 < 4/3). Therefore, Tb is a source for the full state space when 7 > 4/3 and V± are sources for the full state space 
when 7 < 4/3. The numerical experiments have confirmed that Tb is a source in the full state space when 7 > 4/3. 
The numerical simulation is consistent with the decay rates [57] and with the analysis given above that Tb is a source 
for the full state space when 7 = 4/3. 



2. Separable invariant sets 

We set £_ = £ x = 0. We can then show that Ei 1 = E{x)F{t) (i.e., Ei 1 is separable). It follows that Nl + N 2 . = 
F 2 (t) and (q + 3S+ - 2) = F/F. Therefore, setting £_ = S x = B = and taking Q.(t), Q„(t), v(t), F(t) as 
independent functions of r, we obtain 

JV_ = (242) 



F(r)e 91 ( x ) 

E i = * (244) 



and 



Qu = ": iT)n [ T l r , (245) 



l + (7-l)u 2 (r) 

2 



E+(r) = \ (1 - ^ 2 (r) - fi(r) - 0,(r) , (246) 



where 52 ( x ) satisfies the following ODE: 

dg 2 (x) 



dx 

Hence fi(r), n u (r), u(r) and F(t) satisfy the following ODE 

dF(r 



2^3(1 + e 292{x) )e- 9l{x) e- 92{x) . (248) 



dT 

dQ(r) _ 2 7 fi(r) (G + ^ ^ 1, , (1 , 2 



(q + 3E+)F(t) (249) 
(g + l)-i(l-3S + )(l + t; 2 )-l) (250) 

((2- 7 )(l-3E+)t;-2(7-l)i;) (251) 



dr G+ V 7 2 

cfo(r) (1 - « 2 ) 



dr G 

^f 1 = ^2 ((^ + 6S + - 1)« 2 + 3( 9 + 2S + - 1)) (252) 

Therefore, the PDE (211)-(225) have separable solutions (i.e., there is a subset of solutions of the PDE which are 
separable). Because there is no term containing dxEi 1 in these equations, there is no differential equation for g\{x). 
The only equation for g 2 {x) and gi(x) is equation (248). Hence we can take any gi(x), and then determine g 2 (x) 
from equation (248), provided g 2 {x) is a real function for a given real function gi(x). We can choose coordinates to 
set E(x) = 1, which leads to e 91 ^ = vT + e 292 ^ . Thus Eq. (248) becomes 



(1 -J- p^92.\ x )\ I 

2V3 \ =L e' 32 ^ = -2V3a/(1 + e 292 ^)) e 

./1 J- B 202 x) V v ' 



92(x) 



dx Vl + e 292 ( x ) 
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which has real solutions given by 

g 2 (x) = - sinh(2\/3x + C). 

Therefore, it is clear that the spatial dependence is gauge, and that these models are in fact spatially homogeneous. 
Consequently, we can obtain non-trivial spatially homogeneous models (that are difficult to describe due to gauge 
difficulties). Other separable cases can be investigated. 

3. Remarks 

All models have an initial singularity as r — > — oo. In addition, we find that {CI, £_ , 7V X , S x , N- , r} — > as t —* — oo 
for all 7 > 1. In the case 7 > 4/3, the dynamical and numerical analysis indicates that {y, Q u , Q u , £+} — > (and 
fib — > 1) for aZi initial conditions. In the case of radiation (7 = 4/3), the models still isotropize as r — > —00, albeit 
slowly. For 7 < 4/3, {«, f2„, Q„, £+} tend to constant but non-zero values as r — > —00. 

In the invariant set 1; = 0, all models isotropize to the past (for 7 > 1). Thus in the spatially homogeneous case 
with no tilt (with Ei 1 = 0), it follows that there exists an isotropic singularity in all orthogonal Bianchi brane- world 
models in which Clb dominates as the initial singularity is approached into the past (consistent with the results of 
[51]). In particular, Tb is a local source and in general the initial singularity is isotropic. The linearized solution 
representing a general solution in the neighbourhood of the initial singularity in a class of G2 models was given in 
[50] ; it was found that Tb is a local source or past-attractor in this family of spatially inhomogeneous cosmological 
models for 7 > 1. Let us now assume that w ^ (the general case). For the Gi models in the timelike area gauge 
we recover the orthogonally transitive tilting Bianchi type VIo and VIIo models in the spatially homogeneous limit 
with one tilt variable (note that from above other Bianchi models can occur in certain limits). There are no sources 
in the Bianchi type VIo and VIIo models with tilted perfect fluid in GR; the past attractor is an infinite sequence 
of orbits between Kasner points. The Bianchi models have bifurcations at 7 — 2/3 and 4/3 (e.g., the dimension of 
the unstable and stable manifolds of the equilibrium points change) [106]. We note that this is consistent with our 
bifurcation 7 = 4/3 (with 7 — > 27 in brane- world models). 
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VIII. COVARIANT AND GAUGE- INVARIANT LINEAR PERTURBATIONS 



We have found that, unlike standard GR where the generic cosmological singularity is anisotropic, the past attractor 
in the brane for homogeneous anisotropic models and a class of inhomogeous models is a simple BRW model Tb- It 
has been conjectured [50,51] that an isotropic singularity could be a generic feature of brane cosmological models. In 
[73] it was shown that this conjecture is true, within a perturbative approach and in the large-scale and high energy 
regime, through a detailed analysis of generic linear inhomogeneous and anisotropic perturbations [72,76,99] of this 
past attractor Tb by using the full set of linear 1+3 covariant propagation and constraint equations; these equations 
are split into their scalar, vector and tensor contributions which govern the complete perturbation dynamics. 

From a dynamical systems point of view the past attractor Tb for spatially homogeneous cosmological brane- 
world models is a fixed point in the phase space of these models. This phase space may be thought of as an invariant 
submanifold within a higher dimensional phase space for more general inhomogeneous models. The linear perturbation 
analysis can be thought of as an exploration of the neighbourhood of Tb out of the invariant spatially homogeneous 
submanifold. The analysis is restricted to large- scales, at a time when physical scales of perturbations are much 
larger than the Hubble radius, A ^> H . This is not very restrictive, since this is the case for the non - inflationary 
perfect fluid models (that are relevant here) in which any wavelength A < at a given time becomes much larger 
than H^ 1 at early enough times. Hence the analysis is valid for any A as t — > 0. 

Let us restrict the analysis to the case of vanishing background dark radiation {LA — 0) [95] . We define dimensionlcss 
expansion normalised quantities for the shear a a b, the vorticity w a , the electric E a b and magnetic H a b parts of the 
Weyl tensor as follows [99]: 

v _ aab w — L ° a p — ^ ab n-i — ab fo^QA 

^ab — —j-j- , W a — j &ab — ~Jj2 ' nab ~ JJ2 ' 

It is convenient to characterise the vorticity of the fluid flow by using the dimensionless variable 

W* a = acoa (254) 

The appropriate dimensionless density and expansion gradients which describe the scalar and vector parts of density 
perturbations are given by (see [161,162] for details of definitions) 

A a = -D aP , Z* a = ^B a H, (255) 
P ti 

and for the brane -world contributions we define the following dimensionless gradients describing inhomogeneity in 
the non - local quantities 

K = ^U a , Q* a = ^Q a , (256) 

where U a and Q a are defined in equation (27) of [102]. It is also useful to define a set of auxiliary variables corresponding 
to the curls of the standard quantities defined above: 

W* = jjcmiw: , tab = ^curlS^ , (257) 

1 - 1 

Sab = — curlf ah , Hab = JjU ab , (258) 

Q* a = -^curlO: . (259) 

Finally, we use the dimensionless time derivative d/dr = d/d (In a) (denoted by a prime) to analyse the past evolution 
of the perturbation dynamics. 

We use the harmonics defined in [72,76] to expand the above tensors X a ,X a b in terms of scalar (S), vector (V) 
and tensor (T) harmonics Q. This yields a covariant and gauge invariant splitting of the evolution and constraint 
equations for the above quantities into three sets of evolution and constraint equations for scalar, vector and tensor 
modes respectively: 

X = X S Q S (260) 
X a = k- 1 X s Q s a +X v Q v a (261) 
X ab = k- 2 X s Q s ab + k- 1 X v Q v ab + X T Q T ab . (262) 
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The evolution equations for scalar perturbations (suppressing the label S) are given by: 



S' = (37 - 2)S - £ , (263) 

£' = (67 - 3)£ - 37S , (264) 

A' = (37 - 3)A - 7Z* , (265) 

Z*' = (37 - 2)Z* - 3(37 + 1)A - U* , (266) 

Q*' = (3 7 -3)Q* - \U* -67A , (267) 

U*' = (67 - 4)17* , (268) 

and the scalar constraint equations are: 

2£ = 2Z*-3Q*, (269) 

2£ = 6A- 3Q* + f7* . (270) 



The vector and tensor perturbations on large - scales, together with their curls and constraints, are given in [73]. 

Solutions to these three sets of equations are presented in Table I in [73] . This work was extended in [95] by providing 
a complete analysis of scalar, vector and tensor perturbations for different stages of the brane-world evolution. We 
conclude that T\, is stable in the past (as t — > or t — > —00) to generic inhomogeneous and anisotropic perturbations 
provided the matter is described by a non- inflationary perfect fluid with 7- law equation of state parameter satisfying 
7 > 4/3 (and we use the large-scale approximation). In particular, unlike in GR, both the shear and density gradient 
tend to zero at early times if 7 > 4/3. Thus the high -energy models isotropize in the past for realistic equations of 
state when we include generic linear inhomogeneous and anisotropic perturbations. 

We have discussed here only the case of vanishing background Weyl energy density, 14 = 0. This assumption 
considerably simplifies the analysis, but the results remain true for U 7^ 0. Indeed, when W/Ofj still remains 
a past attractor. In other words, our analysis is restricted to the invariant submanifold U = of the larger phase 
space with W/0, but this submanifold is asymptotically stable against hi ^ perturbations. This is consistent with 
the dynamical analysis of the class of spatially inhomogeneous G2 cosmological models in the brane-world scenario 
presented earlier. 
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IX. CONCLUSION 



The governing field equations induced on the brane, using the Gauss-Codazzi equations, matching conditions and 
Z-i symmetry were given [161,162,208,212]. We are particularly interested in the physical case of a perfect fluid with 
p = (7 — l)p satisfying 7 > 4/3, where 7 = 4/3 corresponds to radiation and 7 = 2 corresponds to a massless scalar 
field close to the initial singularity. The dynamical equations on the 3-brane differ from the GR equations in that 
there are nonlocal effects from the free gravitational field in the bulk, transmitted via the projection of the bulk 
Weyl tensor, and the local quadratic energy-momentum corrections, which are significant at very high energies and 
particularly close to the initial singularity. As a consequence of the form of the bulk energy-momentum tensor and 
of Z-i symmetry, it follows [208,212] that the brane energy-momentum tensor separately satisfies the conservation 
equations; i.e., V V T )1V = 0. Consequently, the Bianchi identities on the brane imply that the projected Weyl tensor 
obeys the non-local constraint W^E^ = k^W^S^. The evolution of the anisotropic stress part is not determined 
on the brane; the correction terms must be consistently derived from the higher-dimensional equations. We have 
assumed here that V^ v = 0. 

The asymptotic dynamical evolution of spatially homogeneous brane-world cosmological models close to the initial 
singularity, where the energy density of the matter is larger than the brane tension and the behaviour deviates 
significantly from the classical GR case, has been studied in detail. It has been found that for perfect fluid models 
with a linear barotropic 7-law equation of state an isotropic singularity is a past-attractor in all orthogonal Bianchi 
models and is also a local past-attractor in a class of inhomogeneous brane-world models (for 7 > §)■ The dynamics of 
the class of inhomogeneous brane-world models considered above, together with the BKL conjectures [18-20], indicate 
that a wide class of inhomogeneous brane-world models will have an isotropic initial singularity. 

More detailed information has been obtained from a study of the dynamics of a class of spatially inhomogeneous Gi 
cosmological models with one spatial degree of freedom in the brane-world scenario near the cosmological initial singu- 
larity [57]. Since the normalized frame variables were found to vanish asymptotically, the singularity is characterized 
by the fact that spatial derivatives are dynamically negligible. The local dynamical behaviour of this class of spatially 
inhomogeneous models close to the singularity was then studied numerically [57]. It was found that the area expansion 
rate increases without bound (and hence H — > 00) asr^ —00, so that there always exists an initial singularity. For 
7 > 4/3, the numerics indicate that as t — ► —00 the singularity is isotropic for all initial conditions (indicating that 
Tb is a global past-attractor). In the case of radiation (7 = 4/3), the models were still found to isotropize as r — > — 00, 
albeit slowly. From the numerical analysis we conclude that there is an initial isotropic singularity in all of these G2 
spatially inhomogeneous brane cosmologies for a range of parameter values which include the physically important 
cases of radiation and a scalar field source. The numerical results are supported by a qualitative dynamical analysis 
and a detailed calculation of the past asymptotic decay rates [57]. 

From numerical and dynamical considerations we concluded that the Hubble expansion rate increases without 
bound as logarithmic time r — > —00, and hence there always exists an initial singularity in G 2 cosmologies. In 
addition, the normalized frame variable vanishes asr-^ —00. These conditions allow us to calculate the exponential 
decay rates close to the initial singularity. The constrained evolution system of equations for general inhomogeneous 
(denoted Go) brane-world models was given in [154] in terms of the variables fit,X (in the separable volume gauge 
using Hubble- normalized quantities), where X represents the independent variables corresponding to the shear and 
curvature and additional matter terms, and f2f, = pi,/H 2 , where essentially p^ ~ p 2 . The exponential decay rates in 
the case 7 > 4/3 were calculated [57], and it was found that as r — > —00 the variables (f2ft — 1), X have decay rates that 
depend linearly on {e( 37_1 ) T , e 37r , e 3 ( 7_1 ) T , e 2 ( 37 ~ 4 ) T , e 2r , e( 37 ~ 4 ) T , e 2 ( 37_2 ) T }. There is also evidence of isotropization 
(albeit slowly) in the degenerate case 7 = 4/3. This supports the conjecture that, in general, brane-world cosmologies 
have an isotropic singularity into the past for 7 > 4/3. 

The evolution of models with an isotropic cosmological initial singularity is approximated by the flat self-similar, 
spatially homogeneous and isotropic brane-world BRW model [27,28,83], corresponding to the 'equilibrium state' JF b , 
characterized by fij, = 1, X = 0. It follows immediately that the total energy density p ^ 00 as t ^ - 00 [50], 
so that hi, ~ p 2 dominates as r — > —00, and that all of the other contributions to the brane energy density are 
negligible dynamically as the singularity is approached. Hence close to the singularity the matter contribution is 
given by p tot — -^p 2 = p, b ; p tot — j^(p 2 + 2pp) = (27 — l)p tot , so that the effective equation of state at high densities 
is (27 — 1), which is greater than unity in cases of physical interest. For clarification, let us describe the dynamics 
heuristically. As in GR, the anisotropies (for example) grow into the past as the initial singularity is approached. The 
energy density also diverges in GR, but at a slower rate than the anisotropies and consequently the initial singularity is 
generically anisotropic. However, in brane-world models, at higher energies close to the initial singularity the dominant 
term in the total energy density is effectively p 2 , which now diverges much more rapidly than the anisotropies, so 
that the ratio (i.e., the normalised anisotropies) actually decay, with the result that the anisotropies (and all other 
contributions to the total energy density) are negligible dynamically as the initial singularity is approached. 
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These results are further supported by an analysis of linear perturbations of the BRW model using the covariant and 
gauge invariant approach [73,95]. In particular, a detailed analysis of generic linear inhomogcncous and anisotropic 
perturbations [72,76] of the past attractor T\, was carried out by deriving a full set of linear 1+3 covariant propagation 
and constraint equations for this background. Solutions to the set of perturbation equations were presented [73], and 
it was concluded that Th is stable in the past to generic inhomogeneous and anisotropic perturbations for physically 
relevant values of 7. In addition, it follows immediately that the expansion normalised shear vanishes as the initial 
singularity is approached, and isotropization occurs. 

We have argued that, unlike the situation in GR, it is plausible that typically the initial singularity is isotropic 
in brane- world cosmology. Such a 'quiescent' cosmology [14], in which the universe began in a highly regular state 
but subsequently evolved towards irregularity, might offer an explanation of why our Universe might have began its 
evolution in such a smooth manner and may provide a realisation of Penrose's ideas on gravitational entropy and 
the second law of thermodynamics in cosmology [189]. More importantly, it is therefore possible that a quiescent 
cosmological period occuring in brane cosmology provides a physical scenario in which the universe starts off smooth 
and that naturally gives rise to the conditions for inflation to subsequently take place. Cosmological observations 
indicate that we live in a Universe which is remarkably uniform on very large scales. However, the spatial homogeneity 
and isotropy of the Universe is difficult to explain within the standard GR framework since, in the presence of matter, 
the class of solutions to the Einstein equations which isotropize is essentially a set of measure zero. In the inflationary 
scenario, we live in an isotropic region of a potentially highly irregular universe as the result of an accelerated expansion 
phase in the early universe thereby solving many of the problems of cosmology. Thus this scenario can successfully 
generate a spatially homogeneous and isotropic universe from initial conditions which, in the absence of inflation, 
would have resulted in a universe far removed from the one we live in today. However, still only a restricted set of 
initial data will lead to smooth enough conditions for the onset of inflation. 

Let us discuss this in a little more detail. Although inflation gives a natural solution to the horizon problem of 
the big-bang universe, it requires homogeneous initial conditions over the super-horizon scale; i.e., in general inflation 
itself requires certain improbable initial conditions. When inflation begins to act, the Universe must already be 
smooth on a scale of at least 10 5 times the Planck scale. Therefore, we cannot say that it is a solution of the horizon 
problem, although it reduces the problem by many orders of magnitude. Many people have investigated how initial 
inhomogeneity affects the onset of inflation, and it has been found that including spatial inhomogencities accentuates 
the difference between models like new inflation and those like chaotic inflation. Goldwirth and Piran [96-98], who 
solved the full Einstein equations for a spherically symmetric spacetime, found that small-field inflation models of 
the type of new inflation are so sensitive to initial inhomogeneity that they require homogeneity over a region of 
several horizon sizes. Large-field inflation models such as chaotic inflation are not as severely affected by initial 
inhomogeneity but require a sufficiently high average value of the scalar field over a region of several horizon sizes 
[2,32,138]. Therefore, inhomogeneities further reduce the measure of initial conditions yielding new inflation, whereas 
the inhomogencities have sufficient time to redshift in chaotic inflation, letting the zero mode of the field eventually 
drive successful inflation. In conclusion, although inflation is a possible causal mechanism for homogenization and 
isotropization, there is a fundemental problem in that the initial conditions must be sufficiently smooth in order for 
inflation to subsequently take place [51]. We have found that an isotropic singularity in brane- world cosmology might 
provide for the necessary sufficiently smooth initial conditions to remedy this problem. 

The earliest investigations of the initial singularity used only isotropic fluids as a source of matter. However, the 
study of the behaviour of spatially homogeneous brane-worlds close to the initial singularity in the presence of both 
local and nonlocal stresses indicates that for physically relevant values of the equation of state parameter there exist 
two local past attractors for these brane-worlds, one isotropic and one anisotropic (although the anisotropic models 
are likely unphysical). In particular, Barrow and Hervik [15] studied a class of Bianchi type I brane- world models 
with a pure magnetic field and a perfect fluid with a linear barotropic 7-law equation of state. They found that when 
7 > |, the equilibrium point T\> is again a local source (past-attractor), but that there exists a second equilibrium 
point, denoted PHi, which corresponds to a new brane- world solution with a non-trivial magnetic field, which is also 
a local source. When 7 < |, PH\ is the only local source (the equilibrium point T\> is unstable to magnetic field 
perturbations and hence is a saddle). This was generalized in [54], in which it was shown that for a class of spatially 
homogeneous brane-worlds with anisotropic stresses, both local and nonlocal, the brane-worlds could have either an 
isotropic singularity or an anisotropic singularity for 7 > 4/3. Finally, we note that Bianchi type I and IX models have 
also been studied in Hofava-Witten cosmology in which the fifth coordinate is a S 1 /Z2 orbifold while the remaining 
six dimensions have already been compactified on a Calabi-Yau space, and it was argued [67] that there is no chaotic 
behaviour in such Bianchi IX Hofava-Witten cosmologies (also see [68]). 

It would be of interest to study general inhomogeneous (Go) brane-world models further. The exponential decay 
rates in the case 7 > 4/3 were calculated in [57]. The decay rates are essentially the same as in the G2 case presented 
earlier. This supports the possibility that in general brane-world cosmologies have an isotropic singularity. We hope 
to further study Go brane-world models numerically in the future. The dynamics of inhomogcncous brane-world 
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models, together with the BKL conjectures [18-20], perhaps indicate that a wide class of inhomogeneous brane-world 
models will have an isotropic initial singularity. 

The results might also be applicable in a number of more general situations. For example, in theories with field 
equations with higher-order curvature corrections (e.g., the 4D brane-world in the case of a Gauss-Bonnet term in the 
bulk spacetime [167,187]), the results concerning stability are expected to be essentially valid in regimes in which the 
curvature is negligible (e.g., close to the initial singularity). The single-brane cosmological model can be generalized 
to include stresses other than A in the bulk. The simplest example arises from considering a charged bulk black hole, 
leading to the Reissner-Nordstrom AdS5 bulk metric [11] . The effect of the black hole charge on the brane arises via the 
junction conditions and leads to a modified Friedmann equation, which can cause a bounce at high energies. Another 
example is provided by the Vaidya-AdSs metric, in which a RW brane moves in a 4-isotropic bulk (which is not 
static), with either a radiating bulk black hole or a radiating brane [43,144,147,148]. The metric satisfies the 5D field 
equations with a null-radiation energy-momentum tensor. The modified Friedmann equation has the standard form, 
but with a dark radiation term that no longer behaves strictly like radiation. A more complicated bulk metric arises 
when there is a self-interacting scalar field $ in the bulk. The simplest case, when there is no coupling between the 
bulk field and brane matter, was studied in [10,168,175]. When there is coupling between brane matter and the bulk 
scalar field, then the Friedmann and conservation equations are more complicated. The scalar field could represent 
a bulk dilaton of the gravitational sector, or a modulus field encoding the dynamical influence on the effective 5D 
theory of an extra dimension other than the large fifth dimension [31,35,121,122,134,188,198]. For two-brane models, 
the brane separation introduces a new scalar degree of freedom, the radion. For general potentials of the scalar field 
which provide radion stabilization, 4D Einstein gravity is recovered at low energies on either brane [149,182,183,216]. 
Inflation in brane-world cosmology and the issue of stabilization in models with scalar fields in the bulk has been 
studied [86]. 
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